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AERONAUTIC SYMBOLS 



1. FUNDAMENTAL AND DERIVED UNITS 



Symbol 

Metric 

English 

Unit 

Abbrevia- 

tion 

Unit 

Abbreviation 

T if>n gth 

l 

meter 

m 

foot (or mile) _ 

ft (or mi) 

Time 

t 

second. __ _ 

s 

second (or hour) 

sec (or hr) 

Force 

F 

weight of 1 kilogram 

kg 

weight of 1 pound 

lb 

Pnwpr 

P 

horsepower (metric) _ 


horsepower. _ 

hp 


V 

/kilometers per hour 

kph 

miles per hour. _ _ 

mph 

Speed 

\ meters per second. _ _ 

mps 

feet per second _ 

fps 


w 

9 

m 


M 

S 

S w 

G 

b 

c 

A 

V 

9 

L 

D 

Do 


D t 


D v 

C 


2. GENERAL SYMBOLS 


Weight = mg 

Standard acceleration of gravity =9.80665 m/s 2 
or 32.1740 ft/sec 2 
W 

Mass=— 

9 

Moment of inertia =mk 2 . (Indicate axis of 
radius of gyration k by proper subscript.) 
Coefficient of viscosity 


v Kinematic viscosity 

p Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m“ 4 -s 2 at 15° C 
and 760 mm; or 0.002378 lb-ft“ 4 sec 2 
Specific weight of “standard” air, 1.2255 kg/m 3 or 
0.07651 lb/cu ft 


3. AERODYNAMIC SYMBOLS 


Area 

Area of wing 
Gap 
Span 
Chord 


b 2 


Aspect ratio, -g 
True air speed 
Dynamic pressure, ^ P^ 2 

Lift, absolute coefficient C L =^ 




D 

Drag, absolute coefficient 
Profile drag, absolute coefficient 
Induced drag, absolute coefficient C D =^ 
Parasite drag, absolute coefficient 


C 


Cross-wind force, absolute coefficient 


i w Angle of setting of wings (relative to thrust line) 

11 Angle of stabilizer setting (relative to thrust 

line) 

Q Resultant moment 

12 Resultant angular velocity 

R Reynolds number, p — where Zis a linear dimen- 
sion (e.g., for an airfoil of 1.0 ft chord, 100 
mph, standard pressure at 15° C, the corre- 
sponding Reynolds number is 935,400; or for 
an airfoil of 1.0 m chord, 100 mps, the corre- 
sponding Reynolds number is 6,865,000) 
a Angle of attack 

e Angle of downwash 

oo Angle of attack, infinite aspect ratio 

a< Angle of attack, induced 

<x a Angle of attack, absolute (measured from zero- 

lift position) 

7 Flight-path angle 
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By G. B. Schubauer and H. K. Skramstad 


SUMMARY 

This is an account of an investigation in which oscillations 
were discovered in the laminar boundary layer along aflat plate. 
These oscillations were found during the course of an experiment 
in which transition from laminar to turbulent flow was being 
studied on the plate as the turbulence in the wind stream was 
being reduced to unusually low values by means of damping 
screens. The first part of the paper deals with experimental 
methods and apparatus , measurements of turbulence and 
sound , and studies of transition. A description is then given 
of the manner in which oscillations were discovered and how 
they were found to be related to transition , and then how con- 
trolled oscillations were produced and studied in detail. The 
oscillations are shown to be the velocity variations accompany- 
ing a wave motion in the boundary layer , this wave motion having 
all the characteristics predicted by a stability theory based on 
the exponential growth of small disturbances. A review of this 
theory is given. The work is thus experimental confirmation 
of a mathematical theory of stability which had been in the 
process of development for a period of approximately JfO years , 
mainly by German investigators. 

I. INTRODUCTION 

From previous investigations, both theoretical and experi- 
mental, much is known about the laminar boundary layer. 
Its mean velocity distribution, thickness, and separation 
point may be regarded as so well understood that little 
further study is required. When fluctuations occur in the 
laminar boundary layer, their origin, character, and effect 
on transition from laminar to turbulent flow still present 
formidable problems. It has been observed, for example, 
that the flow in a laminar boundary layer is not steady when 
turbulence is present in the surrounding stream; but, beyond 
the simple observation that velocity fluctuations exist, little 
is known experimentally concerning the behavior of the 
fluctuations. Transition is known to depend on stream tur- 
bulence and fluctuations in the layer are believed to be the 
primary cause of transition, but just how the three are inter- 
related has never been completely understood. 

Experimental evidence indicates, in some cases at least, 
that transition results from separation of the boundary layer 
caused by adverse pressure gradients; and the view has been 
taken by Taylor (reference 1) and other investigators that 
the local pressure gradients accompanying velocity fluctua- 
tions bring about transition in this way. Isotropic turbu- 
lence in the surrounding flow involves velocity fluctuations 


for which the pressure gradients are known when the in- 
tensity and the scale of the turbulence are specified, and 
relations involving this pressure gradient and transition 
Reynolds number that are in fair agreement with experiment 
have been proposed. These relations are restricted to a par- 
ticular type of disturbance but are useful for expressing 
experimental results on the effect of isotropic turbulence. 

On the purely theoretical side, attempts have been made 
to account for transition by stability considerations in 
which the growth or decay of a small disturbance is deter- 
mined by the kinematic viscosity of the fluid, the thickness 
of the boundary layer, the speed outside the layer, and the 
frequency or the wave length of the disturbance. Important 
contributions to stability theory have been made by Tollmien 
(references 2 and 3) and Schlichting (references 4 to 7). 
Experimental evidence has so far yielded only scant support 
of the stability theory; and this support has come from 
experiments in which the flow was known to be of an unstable 
type, namely, the flow in a divergent channel (reference 8), 
in a wake at low speeds (reference 9), and in an acoustically 
sensitive jet (reference 10). No support has come from past 
wind-tunnel experiments, in which the connection between 
the disturbance and transition appears to be mainly through 
the magnitude of the disturbance. In short, the accumu- 
lated wind-tunnel results have led to the general opinion 
among experimentalists that the pressure gradients accom- 
panying any disturbance in the surrounding flow cause tran- 
sition when such gradients become sufficiently large. In 
contrast, stability theory requires that the boundary layer 
be either stable or unstable, depending merely on the fre- 
quency or the wave length of the disturbances present; 
hence, the nature of the disturbance rather than its magni- 
tude is the essential quality. Experiment and theory have 
therefore led to divergent points of view and neither has 
fully explained the nature of transition. 

In the present experiment it is hoped that the two points 
of view have been brought into closer accord for one of the 
simplest, but practical, boundary-layer problems — that of 
the Blasius distribution near the surface of a flat plate. The 
experiment was originally intended to be a general study of 
transition on a flat plate at zero pressure gradient under 
conditions of low stream turbulence, one purpose being to 
trace the course of transition as near as possible to vanishing 
turbulence. This purpose was carried out, but the appear- 
ance at low turbulence of an oscillation in the laminar 
boundary layer opened a new field of investigation and led 
to an extensive study of the phenomena termed “laminar- 
boundary-layer oscillations.” Investigations revealed that 
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the oscillations were the result of amplification of small 
disturbances in accordance with stability theory. When 
the oscillations were artificially produced, their character- 
istics could be studied in detail and compared with the 
characteristics prescribed by theory. 

It seems evident from the results of this investigation 
that transition involves two things: One is the disturbance 
itself, its origin, and its behavior in the boundary layer; the 
other is the effect of the disturbance, either produced in the 
boundary layer or entering from without, in bringing about 
turbulent flow. Sufficiently small disturbances cannot pro- 
duce transition, but a small disturbance may grow according 
to stability theory until it is sufficiently large to cause the 
flow to become turbulent. The growing process is the link 
between the stability theory and transition. 

In most of the earlier investigations, observations of the 
transition point were made but little or no effort was made 
to observe disturbances in the boundary layer before tran- 
sition. The first known observations of velocity fluctuations 
in the laminar boundary layer were made by Dryden (refer- 
ence 11) but, in his experiment as in most others, the initial 
disturbances were so large that transition occurred after only 
a small growth. Nikuradse (reference 12) attempted to test 
the stability theory by producing sinusoidal fluctuations in 
the boundary layer near the leading edge of a flat plate in 
water; but the results were inconclusive, partly because the 
artificial fluctuations, as well as general disturbances, were 
too large and partly because observations were limited to 
the determination of transition points. 

The results obtained in the present investigation can be 
attributed largely to the low level of turbulence in the wind 
stream and to the careful study made of velocity fluctuations 
throughout the laminar boundary layer. Hot-wire ane- 
mometers of high sensitivity and of special design made such 
an investigation possible. 

The work was conducted at the National Bureau of Stand- 
ards under the sponsorship and with the financial assistance 
of the National Advisory Committee for Aeronautics. The 
investigation of transition on a flat plate at low turbulence 
was suggested by Dr. Hugh L. Dryden, and the authors 
wish to acknowledge his assistance and many valuable sug- 
gestions in connection with the experimental program. 


II. SYMBOLS 


X 

y 

u 0 

u 

u 

v 

V) 


u ' , v ' , w' 

Ct 

/3 r =2-7r/ 

/ 

a= 2 t/\ 
X 


P 

M 

v 

<Z 

8 

5* 


R 

Rr 


A 


distance from leading edge of flat plate 
distance from surface of flat plate 
mean velocity outside boundary layer 
mean velocity at a point in boundary layer 
instantaneous ^-component of fluctuation velocity 
instantaneous y-component of fluctuation velocity 
instantaneous fluctuation velocity perpendicular to 
U 0 and parallel to surface of flat plate 
root-mean-square values of u , y, and w, respectively 
wave velocity 

oscillation frequency 


wave length 
amplification coefficient 
density 
viscosity 

kinematic viscosity (m/p) 
dynamic pressure 
boundary-layer thickness 
boundary-layer displacement thickness 


5*= 1 .72 ^ ^ j for Blasius velocity distribution 

5* = 0.341 8 relation used by Tollmien and 
Schlichting (references 2, 4, and 5) 
Reynolds number ( U 0 8*/v ) 

^-Reynolds number ( U 0 x/v ) 

7?=1.72 a /R^ for Blasius velocity distribution 

/ d U 8~\ 

Karman-Polilhausen parameter ( -Vy - j 


III. APPARATUS AND METHODS 

1. WIND TUNNEL 

The present investigation was conducted in the 4K-foot 
wind tunnel at the National Bureau of Standards. The 
general layout of the tunnel is shown in figure 1. The flat 
plate was located vertically in the test section with the leading 
edge 6 feet from the upstream end. In order to reduce vibra- 
tion, the test section of this tunnel is supported directly from 
the foundation and is joined to the rest of the circuit by only 



Figure 1 .— Elevation view of 4K>-foot wind tunnel. 
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a sponge-rubber seal at each end. The area reduction from 
settling chamber to test section is 7.1:1. The guide vanes 
ahead of the settling chamber were made finer than the others 
in order to reduce the scale of the turbulence and to permit 
as much reduction in turbulence through decay as possible. 
Since directional fluctuations in the horizontal plane were 
found to be large, closely spaced straighteners were placed 
at right angles to and on top of the fine guide vanes. This 
combination was in effect a honeycomb and resulted in great 
improvement in the steadiness of the stream. Further reduc- 
tion of turbulence was obtained by installing damping screens 
in the settling chamber. Without screens, the turbulence in 
the test section was 0.27 percent at 100 feet per second and, 
with seven screens, was 0.032 percent. 

2. FLAT PLATE 

The flat plate was a commercial sheet of aluminum % inch 
thick, 4 }& feet wide, and 12 feet long. The leading edge, 
which was symmetrical and pointed, was formed by the 
intersection of circular arcs tangent to the original surface 
4 inches from the leading edge. The surface was left in its 
original condition and had a mirrorlike finish, marred slightly 
by small scratches barely perceptible to the fingertips. Wav- 
iness was quite perceptible when viewed near grazing inci- 
dence. A surface gage showed variations from approxi- 
mately 0.01 to 0.02 inch over distances from 1 to 2 feet. 
These were bends in the plate evidently produced by the 
rolling process. 

The plate was bolted along the top and the bottom edges 
to one side of 3-inch steel channels, which were in turn bolted 
to the floor and to the ceiling of the tunnel. The flanges 
away from the plate served as rails at the floor and at the 
ceiling for a carriage on which all exploring apparatus was 
mounted. The carriage was simply a steel plate % inch thick 
and 6 inches wide, with runners and guides at the ends in 
contact with the rails. This plate was thus parallel to the 
wind and 3 inches from the flat plate. All explorations in 
the boundary layer were made 25 inches or more ahead of 
the carriage to avoid the local pressure field. The carriage 
was propelled by hand from the outside of the tunnel by a 
sprocket and chain. 

3. CONTROL OF PRESSURE DISTRIBUTION 

In order to control the pressure gradient along the flat 
plate, the cross section of the working chamber was varied by 
adjustable auxiliary walls on the vertical sides of the cham- 
ber. These walls were aluminum sheets }{ 6 inch thick and 
22 inches wide extending from opposite the trailing edge of 
the plate to a distance 4 feet ahead of the leading edge. 
The sheets, mounted on screws spaced 9 inches and threaded 
through the tunnel walls, were adjustable to and from the 
sides of the tunnel. In order to prevent discontinuity at the 
top and at the bottom edges, the sheets were backed by a 
stretched rubber diaphragm; the combination of aluminum 
sheet and rubber formed a continuous flexible wall that 
could be warped or displaced as desired. The maximum 
range of displacement of these auxiliary walls was about 
4 inches. 


By this means the pressure could be made to rise, to fall, 
or to remain constant along the plate. (See curves A, B, and 
C, fig. 30.) For most of the experiment the walls were set 
for zero pressure gradient. 

The provision for warping the walls into a bulge or a 
hollow did not prove so useful as was anticipated because 
the effect of the warp usually extended over too great a dis- 
tance to give the desired result. When it was desired to 
produce a sharp pressure rise or fall, other devices were used. 
For example, an airfoil extending from floor to ceiling near 
the plate produced a sharp pressure fall followed by an 
abrupt rise. A pressure fall without a rise was easily pro- 
duced by proper blocking of the stream. 

Although measurements were taken on only one side of the 
flat plate, the auxiliary walls on each side of the tunnel were 
given the same shape for symmetry. Some asymmetry 
nevertheless existed because of the presence of the carriage 
and the measuring equipment on only one side of the plate. 
Although this resulted in greater blocking of the stream on 
the working side, it was found that some excess blocking here 
was necessary to direct the flow at a slight angle to the lead- 
ing edge of the plate so that the stagnation point was dis- 
placed slightly from the sharp edge to the working side of 
the plate. The transition point was unaffected by the posi- 
tion of the stagnation point so long as this point was not 
displaced to the opposite side of the plate. The slight excess 
blocking was therefore maintained at all times to prevent 
directional variations accompanying the turbulence in the 
stream from ever carrying the stagnation point to the 
opposite side. 

4. MEASUREMENT OF PRESSURE DISTRIBUTION AND DETERMINATION OF 
TRANSITION POINT 

A small pitot-static head was arranged as shown in figure 
2 with the impact tube in contact with the surface and the 
static tube parallel to the stream % inch from the surface. 
This combination was carried on an arm attached to the car- 
riage 28 inches to the rear and could slide fore and aft with 
forklike guides and impact tube always in contact with the 
surface. Both tubes were made from thin-wall nickel tubing 



Figure 2.— Pitot-static surface tube used to measure pressure distribution and position of 

transition. 
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with 0.04-inch outside diameter. The impact tube was 
flattened on the end to form a slit 0.007 inch wide, the center 
of which was 0.006 inch from the surface when the tube was 
in contact. The static tube had four 0.008-inch holes drilled 
through the wall, 8 diameters from the closed end. 

Only the static tube of this instrument was used when 
pressure distributions were measured and both tubes were 
used when transition points were determined. Traverses 
were always horizontal and usually made along the center 
line of the plate. The static pressure was thus measured 
/4 inch from the surface, which for most of the surface corre- 
sponded to a position outside the boundary layer. A large 
variation of pressure normal to the surface was found near 
the leading edge where the pressure distribution was deter- 
mined by the shape of the leading edge. Also, pressure 
effects caused by waviness of the surface depended on the 
distance from the surface. Since the former variations were 
confined to the first 6 inches from the leading edge and the 
latter were small, variations in pressure in the K-ineh distance 
normal to the surface were neglected. 

The detection of transition by the surface-tube method 
depends on the variation with x of impact pressure near the 
surface, which corresponds closely to the variation of shear- 
ing stress at the surface. This impact pressure decreases 
with distance from the leading edge to the beginning of 
transition, then rises through the transition region, and again 
falls in the turbulent region. In the present work, the point 
of minimum pressure was taken as the beginning of transi- 
tion and the point of maximum pressure following the rise 
was taken as the point where turbulence was fully developed. 

5. DETERMINATION OF BOUNDARY-LAYER THICKNESS 

Velocity distributions across the boundary layer were de- 
termined by traversing normal to the surface with a flat- 
tened impact tube similar to the surface tube. This tube 
was carried on an arm similar to the one shown with the hot- 
wire anemometer attached in figure 3. Both the micrometer 



3" 

Figure 3.— Apparatus for traversing normal to the surface. Hot-wire head attached. 


screw A and the fulcrum B were held in contact with the 
surface of the plate by the rigid arm extending rearward to 
the carriage, and motion of the impact tube to and from the 
surface was obtained by rocking the arm about the fulcrum 
by means of the screw. The velocity distributions were 
used to determine <5*. When the pressure gradient was 


zero, the measured values of <5* were in agreement with those 
calculated by the Blasius formula 

Whenever the value of 5* was needed for presenting results 
obtained at zero pressure gradient, the value calculated by 
this formula was used. 

Tollmien and Schlichting (references 2, 4, and 5) used a 
constant in the above formula equal to 1.73 for computing 
5*. More recent values of the constant are 1.7207 given by 
Dryden (reference 11) and 1.7208 given by Goldstein (refer- 
ence 13). A value of 1.72 is therefore used in the present 
work. 

6. HOT-WIRE TURBULENCE EQUIPMENT 

A variety of apparatus, such as amplifiers, bridges, poten- 
tiometers, oscillators, oscillographs, and numerous types of 
hot-wire anemometer, comes under the heading of hot-wire 
equipment for measuring turbulence. Amplifiers with 
attendant power supplies and liot-wire anemometers are 
usually special equipment, whereas the other equipment is 
of standard commercial design. Amplifiers used in turbu- 
lence work and circuits to compensate for the thermal lag 
of hot wires have, however, been standardized to a certain 
extent and their essential features are described in the 
literature. Two amplifiers, which were designed and 
built by W. C. Mock, Jr., were used in the present investi- 
gation. One of the amplifiers is described in reference 14; 
the other was a newer and more portable type with about 
the same frequency response but using capacitance compen- 
sation instead of inductance compensation. Both types of 
compensation and the requirements to be met by each type 
are described in reference 15. The theory of compensation 
is given in reference 16. When properly compensated, the 
over-all response of wire and amplifiers was uniform from 3 
to about 2000 cycles per second. 

In all cases hot-wire anemometers must be calibrated by 
measuring the voltage across the wire at known airspeeds;* 
or, if the anemometer is sensitive to direction changes, it 
must be calibrated by measuring the voltage at several angles 
to the wind. For this purpose a potentiometer is necessary. 
A Wheatstone bridge must be used to obtain the resistance 
of the wire at air temperature. Throughout the present 
work, both in calibration and in use, the heating current 
was held constant and only the temperature and the volt- 
age across the wire were allowed to vary. This condition 
was maintained during rapid velocity fluctuations by a suffi- 
ciently large choke coil in series with the wire. Amplifiers 
were calibrated by applying to them a known alternating 
voltage from an oscillator. By use of the calibration data, 
velocity fluctuations could be calculated from the amplified 
and properly compensated voltage fluctuations across the 
hot wire. A cathode-ray oscillograph was used for visual 
observation of fluctuations and also for making photographic 
records with a moving-film camera. When observations 
were simultaneously taken on two wires, an electronic switch 
was used in conjunction with the oscillograph. 
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7. HOT-WIRE ANEMOMETERS 

The hot-wire anemometers used here may be divided into 
two classes: (1) Those used in the free stream to measure 
turbulence; (2) those used in the boundary layer to study 
oscillations, transition phenomena, and turbulence. An- 
emometers of class (1) were designed for high sensitivity and 
for freedom from vibrational effects. Anemometers of class 
(2) were designed for working near a surface with a mini- 
mum of interference and with as little vibrational motion 
relative to the surface as possible. Both types used plati- 
num wire 0.00025 inch in diameter, obtained by etching the 
silver coating from Wollaston wire. The etched wire was 
soft-soldered to prongs made either from fine sewing needles 
or from fine copper wire. Platinum wires of this size were 
adequately sensitive in lengths no greater than Ye inch and, 
in the range of operating conditions, had time constants rang- 
ing from 0.0005 to 0.002 second. 

Anemometers of class (1) were made with a single wire 
normal to the wind for measuring u' or with two wires set 
at an angle to the wind in the form of an X for measuring 
v' and w ' . Heads of these two types are shown in figure 4. 



r 


Figure 4.— Hot-wire heads for the measurement of free-stream turbulence. 

Both have platinum wires (not visible in the photograph) 
across the tips of the prongs; the four-prong head lias the two 
wires forming an X. In the X -arrangement the wires lie 
as nearly as possible in a plane without touching and sub- 
tend an angle of 60° with each wire at an angle of about 30° 
to the wind. Braces of silk thread cemented across the 
prongs about % inch from the tips were necessary to prevent 
wire breakage. All wires were about Ye inch long. These 
arrangements of wires take advantage of the directional 
characteristics of a wire in order to obtain sensitivity to 
the component desired. The sensitivity of a wire to direc- 
tion is known to depend on the angle to the wind, being zero 
when normal and again when parallel. Sensitivity to u 
is a maximum when the wire is normal. Since v and w are 
small, their principal effect is to cause fluctuation in direction 
of the stream when added vectorially to U and the effect on 
the magnitude of the instantaneous velocity is insignificant. 
A wire normal to the wind responds thus to ^-fluctuations 
rather than to v or w. 

An X-arrangement of two identical wires, each making the 


same angle to the wind, with voltage opposed produces a 
resultant voltage change only when the wires are differentially 
cooled. This arrangement is then insensitive to u but highly 
sensitive to direction changes in the plane of the wires and 
hence to v or to w depending on the orientation of the plane. 
When fluctuations are large, isolation of single components 
in this way is not complete and voltage changes are not 
exactly proportional to velocity changes. If the fluctuations 
do not exceed 5 percent of the mean speed, errors from these 
sources are believed to be less than 1 percent. 

The mountings for the hot-wire heads are shown in figure 5. 
The central member to which the head is attached is held 



Figure 5.— Rubber-suspension type mountings with hot-wire heads attached. 


within the cylindrical tube by a rubber-band suspension and 
in this way is isolated from vibration from the tunnel walls. 
Vibrational effects were not troublesome after a sufficiently 
“soft” suspension of this sort was used. The mounting for 
the X- wires contains an angle-changing device for calibrating 
the wires. 

In the boundary layer only the ^-component was measured, 
except when one attempt was made to discover whether the 
boundary-layer oscillations had a w-compopent. A study 
of v was considered desirable but no appropriate hot-wire 
head of sufficiently small dimensions was available for the 
purpose. For anemometers of class (2), therefore, most 
attention was given to hot-wire heads with the wire normal 
to the wind and parallel to the surface. The wire length was 
usually 0.04 inch. For traversing along the surface with the 
wire at a fixed distance, the wire was attached at the tips of 
prongs extending about % inch forward from a small celluloid 
sled, which was held against the surface by wire springs from 
a brass tube attached to the carriage. A sled assembly with 
two pairs of prongs for supporting two wires at different 
distances from the surface is shown in figure 6. Another 
sled was also used with two hot wires making a V in a plane 
parallel to the surface. The latter arrangement was sensi- 
tive to w. The sled was believed to be the best possible type 
of hot-wire head to prevent relative motion between the wire 
and the surface. Interference effects were found negligible. 

For traversing normal to the surface, the hot wire was 
supported on the tips of steel sewing needles extending into 
the boundary layer from the arm of the traversing apparatus 
shown in figure 3. This arrangement was made sufficiently 
rigid to be free from vibration troubles but was not appro- 
priate for traversing parallel to the surface. A combination 
of this apparatus with a sled was used when it was desired 
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Figure G. — Sledlike hot-wire head arranged for two hot wires. 


to keep one wire at a fixed distance from the surface while 
moving the other in and out from the surface. 

IV. REDUCTION OF STREAM TURBULENCE AND EFFECT ON 
REYNOLDS NUMBER OF TRANSITION AT ZERO PRESSURE 
GRADIENT 

In the present investigation, effects of high stream turbu- 
lence were of little concern. Emphasis was placed rather on 
reducing the turbulence as far as practicable in order to 
study the boundary layer when it was as little disturbed from 
the stream as possible. The turbulence was successively 
reduced by placing various numbers and combinations of 
damping screens in the settling chamber. For each addi- 
tion of a screen and for every recombination of screens, 
measurements were made of u'/U 0 , v'/U 0} and w'/U 0 in the 
working chamber and the values of R x at the beginning and 
at the end of transition were determined. 

No correction for wire length according to the method given 
in reference 17 was made, since the appropriate scale for v' 
and w' could not be determined. If the probable order of 
magnitude of the scale is considered, the correction is small 
and is believed to be unimportant in view of the accuracy 
attainable in the measurement of low turbulence. 

Screens of fine mesh and small wire were used in order 
that the turbulence produced by the screen itself would 
decay rapidly. Each screen completely covered the cross 
section of the settling chamber and, when more than one 
screen was used, a spacing of 6 inches or more was allowed 
between them. Preliminary surveys with the hot-wire 
anemometer showed that the turbulence was uniform over 
the usable cross section of the stream and that little or no 
decrease in turbulence occurred in the length of the working 
chamber. The turbulence could therefore be specified by 
hot-wire measurements made at one point in the working 
chamber. The measurements showed that v' and w' were 
nearly equal to each other but generally greater than u ' . 
This nonisotropic condition is known to be caused by the 
contraction of the stream from the large area of the settling 
chamber to the smaller area of the working chamber. The 
more the turbulence was reduced by the addition of screens, 
the less this difference became and, finally, near the lowest 


turbulence u' was actually larger than v' and w r . In all 
cases the turbulence increased with the wind speed, an effect 
that could be partly accounted for bv the decreasing resist- 
ance coefficient of the screens with increasing speed. 

TABLE I 

TURBULENCE WITH SIX DAMPING SCREENS IN SETTLING 

CHAMBER 

[Measurements at single station in working chamber] 


77« 

u' 

r' 

w' 

yy(u ,2 +0' 2 +M> /2 ) 

<-'0 

Uo 

Uo 

Uo 

Uo 

(fps) 

(percent) 

(percent) 

(percent) 

(percent) 

30 

0.019 

0.011 

0. 012 

0.0145 

40 

.021 

.020 

.019 

.0200 

50 

.024 

.023 

.020 

.0224 

60 

.027 

. 026 

.021 

.0248 

70 

.030 

.030 

.026 

.0287 

80 

.033 

.035 

.032 

.0334 

90 

.035 

.039 

.037 

. 0371 

100 

.037 

.042 

.041 

.0401 

110 

.040 

.044 

.045 

.0430 

120 

.045 

.046 

.047 

.0459 


As the number of screens was increased, a reduction in 
turbulence was obtained with each additional screen but by 
successively decreasing amounts; and the reduction by this 
means appeared to have reached a practicable limit when six 
screens were installed. The measured values with six 
screens are given in table I. It was believed at first that a 
level had been reached nearly equal to the turbulence pro- 
duced by the screens themselves. A finer screen than any 
of the others was therefore added downstream from the 
six screens already present. Since no screens with finer wires 
than those already used were commercially available, a 
screen made from silk bolting cloth was installed. Some 
further reduction was obtained but not so much as was 
expected. No attempt was made at further reduction. The 

TABLE II 

TURBULENCE WITH SIX DAMPING SCREENS FOLLOWED 
BY ONE BOLTING CLOTH IN SETTLING CHAMBER 


[Measurements at single station in working chamber] 


Uo 

(fps) 

i/' 

Vo 

(percent) 

p' 

Uo 

(percent) 

w' 

Uo 

(percent) 

(u' 2 +P' 2 +tfl' 2 ) 

Uo 

(percent) 

30 

0.018 

0.012 

0. 009 

0. 0135 

40 

.019 

.012 

.010 

.0142 

50 

.021 

.014 

.015 

.0169 

60 

.026 

.018 

.016 

.0204 

70 

.029 

.020 

.017 

.0226 

80 

.033 

.022 

.018 

.0251 

90 

.037 

.026 

.020 

.0286 

100 

.040 

.029 

.023 

.0315 

110 

.044 

.034 

.025 

.0352 


measured values of the turbulence at this lowest level are 
given in table II. The ^-Reynolds numbers of transition 
found during the various stages of reduction are shown in 
figure 7. The highest turbulence shown in the figure (0.342 
percent) was obtained by placing a square-mesh grid of 
K-incli rope, spaced 6 inches, in the settling chamber a few 
inches downstream from the set of six wire screens. The 
transition point found by Hall and Hislop (reference 18) 
for their lowest turbulence is indicated in figure 7. By 
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definition, their transition point corresponds to the begin- 
ning of transition in the present paper. 

Two questions are raised by the foregoing results: (1) 
Why was the addition of the bolting-clotli screen so inef- 
fective in further reducing the turbulence, and (2) why was 
transition on the flat plate unaffected by reduction of tur- 
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Figure 7.— Effect of turbulence ona-Reynolds number of transition. Flat plate; zero pressure 

gradient. 

bulence below about 0.08 percent? The first question is 
answered in the following section. The second question is 
more difficult and an attempt will be made to answer it only 
after boundan^-Iayer oscillations are considered. (See 
section IX.) 

Y. RESIDUAL TURBULENCE AND NOISE 

While the measurements of the lowest turbulence were in 
progress, it was noted by the trace on the screen of a cathode- 
ray oscillograph that the velocitj^ fluctuations did not have 
the random character usually ascribed to turbulence. 
Vibration of the hot-wire mounting was suspected but was 
ruled out since changes in the softness of the mounting 
failed to change the appearance of the trace. Consideration 
was next given to the effect of wind-tunnel noise and led to 
the startling conclusion that a reasonably loud sound could 
involve particle velocities sufficient to produce an apparent 
turbulence of the order of that actually found. For example, 
a calculation based on a plane wave with an intensity of 
105 decibels above a base level of 10 -9 erg per square centi- 
meter per second showed a root-mean-square particle 
velocity of 0.028 foot per second. At a wind velocity of 100 
feet per second, this is equivalent to an apparent turbulence 
of 0.028 percent, while the measured turbulence was 0.0315 
percent. A few measurements of intensity with a sound 
meter connected to a crystal microphone showed that the 
noise level was probably between 105 and 110 decibels at 
100 feet per second. Since 110 decibels is equivalent to an 
apparent turbulence of 0.051 percent, the noise could 
account for the entire hot-wire reading at this speed. At the 
lower speeds, however, the noise intensity appeared insuffi- 
cient to account for the measured turbulence. Since the 

843784 — 50 2 


wave pattern in the tunnel was very complex, calculations 
based on plane waves can be at best only rough approxima- 
tions. Furthermore to avoid the direct effect of the wind on 
the microphone, it was necessary to measure the sound level 
at a distance from the hot wire. 

An analysis was made by a wave analyzer of the spectrum 
of the noise from the tunnel picked up in the tunnel control 
room by the crystal microphone. A similar analysis was 
then made of the output from the hot wire in the wind 
tunnel. In both cases the wind velocity was 80 feet per 
second. A comparison of the two spectrums is shown in 
figure 8. Since the distribution of intensity with frequency 
is compared here, uncertainties in absolute intensity are of 
minor importance. The similarity of the curves indicates 
that, at a wind velocity of 80 feet per second, a large part 
of the hot-wire output is derived from noise. 



Figure 8.— Distribution with frequency of u-component of turbulence and sound. Wind 
speed, 80 feet per second; u'JUo (total) =0.033 percent. 

The values of turbulence in tables I and II probably 
contain a significant contribution from noise, especially at 
the higher velocities. A large part of the increase with 
velocity is probably due to noise. Also, a u' greater than 
v' and w f is evidence that not all the contributing fluctua- 
tions were carried with the stream through the entrance cone. 
Since most of the noise comes from the propeller, it is obvious 
that damping screens can reduce only the true turbulence. 
On the other hand, the addition of screens must increase the 
noise because of the higher propeller speed needed to maintain 
a given wind velocity through the added resistance in the 
screen. These facts afford a satisfactory explanation of the 
effective limit of damping screens in the present tunnel. 

Figure 8 also suggests some information concerning the 
answer to the second question posed in section IV regarding 
the effect of turbulence on transition. Since the spectrum 
of this apparent turbulence is unlike that of isotropic tur- 
bulence (reference 19), such turbulence will probably affect 
transition in a different manner from isotropic turbulence. 
Velocity fluctuations from sound involve pressure changes 
not related to those of turbulence and involve motions cor- 
related over great distances. Since the proportion of noise 
to turbulence was increased with each reduction of tur- 
bulence, anomalous transition effects are to be expected. 
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VI. LAMINAR-BOUNDARY-LAYER OSCILLATIONS, ZERO 
PRESSURE GRADIENT 

1. FIRST EVIDENCE OF OSCILLATIONS 

When the stream turbulence had been reduced to nearly 
its lowest level with the six wire damping screens, it was 
decided to investigate the velocity fluctuations in the 
laminar boundary layer itself. For this purpose one of the 
sledlike heads with a wire sensitive to u was arranged to 
slide along the surface with the wire 0.023 inch from the 
surface. With greatly reduced stream turbulence, it was 
expected that the slow irregular fluctuations, reported for a 
much more turbulent stream in reference 11, would be cor- 
respondingly reduced. Such fluctuations were, in fact, 
almost nonexistent; but, as the wire was moved downstream 
through the boundary layer, a regular oscillation appeared, 
weak at first but with increasing amplitude as the distance 
downstream increased. Just ahead of the transition, bursts 
of very large amplitude occurred and, at the initial point 
of transition, these bursts were accompanied by a breaking 
into an irregular high-frequency fluctuation characteristic 
of turbulence. Records of the oscillations were made by 
photographing the screen of the cathode-ray oscillograph 
with a moving-film camera. Two sets of such records are 
reproduced as figures 9 and 10. Figure 9 shows the pro- 


same phenomenon at a lower wind speed and at a greater 
distance from the leading edge of the plate is shown in 
figure 10. Here the amplifier gain (magnification) was 
adjusted to keep the recorded amplitude the same from posi- 
tion to position in order to bring out the weak as well as the 
strong oscillations. Attention is called to the lower oscilla- 
tion frequency corresponding in this case to the lower speed 
and to the thicker boundary layer. 



■Stf ^\AaA 




Figure 9.— Oscillograms of u-fiuctuations showing laminar boundary-layer oscillations in 
boundary layer of flat plate. Distance from surface, 0.023 inch; Uo=80 feet per second; 
time interval between dots, Ho second. 


gressive development of the oscillations. The traces at 
6.25 and 6.5 feet already show the bursts of large-amplitude 
oscillations interspersed with highly disturbed turbulent 
motion. The final record at 8 feet indicates a completely 
turbulent boundary layer. The frequency of the oscilla- 
tions may be judged by the Ko-second timing dots. The 



Figure 10.— Oscillograms of ^-fluctuations showing laminar boundary-layer oscillations in 
boundary layer of flat plate. Distance from surface, 0.023 inch; Uo=53 feet per second; 
time interval between dots. Ho second. 

It was thought at first that these oscillations might possi- 
bly be due to vibration of the wire that gave rise to relative 
motion between the wire and the surface. This was soon 
ruled out since changes in the hot-wire head and support 
failed to cause any change in the oscillations. Furthermore, 
the oscillations appeared to be definitely connected with 
transition since the zone in which they occurred always pre- 
ceded transition and moved with it fore and aft along the 
plate as the wind speed was varied. It could not be said 
with certainty that the oscillations were the cause of transi- 
tion, since there existed the possibility that the boundary 
layer became shock-excited by transition occurring a short 
distance downstream, giving rise to an oscillation that was 
possibly the result of transition rather than the cause. 
The latter explanation was ruled out by removing transition 
with an abrupt pressure drop and yet leaving the oscilla- 
tions totally unchanged at an upstream position. The alter- 
native conclusion, with no evidence to the contrary, was that 
transition resulted from the growing of the oscillations to 
the point where the boundary layer was so highly disturbed 
that transition occurred. 

It was soon found that the hot wire could be placed at 
any cross section of the boundary layer to pick up the oscilla- 
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tions, although the amplitude quite evidently varied with 
position. The frequency remained the same throughout the 
section except very near the surface where there was a sug- 
gestion of frequency doubling on the low-velocity part of 
the cycle. This effect at the 4-foot position is shown in 
figure 11. Since the hot wire responds to negative velocities 
as though they were positive, the doubling was believed to 
indicate that the amplitude of the oscillation was sufficiently 
high to reverse the direction of flow during half of each cycle. 

A V-wdre, sensitive to revealed a ^-component in the 
oscillations. Little work was done with this wire and the 
relative amplitudes of u and w were not compared. With 
the one exception cited here, all work was done with wire 
sensitive to u. Because of experimental difficulties no at- 
tempt was made to detect the ^-component, although it is 
shown by theoretical considerations in section VI, 2, that 
a r-component must exist. The presence of the ^-component 
thus indicates that the oscillations were three-dimensional. 

Because the effects of noise were so much in evidence, it 
might be supposed that the oscillations were a resonant 
acoustic phenomenon arising from some frequency in the 
sound spectrum. Among the evidence against this supposi- 
tion is that presented in figure 12, which shoves the effect of 
a tenfold increase in turbulence produced by the rope grid 
in the settling chamber downstream from the six damping 

Wire / 

0.004 in. from surface 

.25 ^ ^ 


screens. With the higher turbulence, it is seen that oscilla- 
tions are present where they were absent before. It thus 
appears that disturbances in the stream, possibly acoustic 
as well as turbulent, give rise to oscillations which are not 
themselves sound waves. 

It was believed almost from the start that the oscillations 
were amplified disturbances, their purity resulting from 
selective amplification of a single frequency or, at most, a 
narrow band of frequencies. Theories had been advanced to 
account for just such an amplification but were often ignored 
in experiment because no conclusive evidence in their support 
had been found. In order to test whether these were in 
reality the amplified oscillations predicted by theory, fre- 
quencies were determined from numerous oscillograms, taken 
at as many positions and speeds as possible and plotted in 
an appropriate manner on a diagram derived from the 
theoretical results of Schlichting. The diagram with the 
experimentally determined points is shown in figure 13. 
This figure will be made clear by the theoretical discussion 
in section VI, 2, and will be described in detail in section 
VI, 3. For the present, it is sufficient to observe that the 
experimental points are distributed along branch II of the 
theoretical curve, which is the region in which the points 
should fall if the initial disturbances, out of which oscilla- 
tions grew, have received approximately their maximum am- 
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Figure 11.— Oscillograms of w-fluctuations showing evidence of frequency doubling near surface. Uo= 105 feet per second; time interval between dots, Ho second. 
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Figure 12.— Oscillograms of u -fluctuations showing effect of stream turbulence on boundary-layer oscillations. Distance from surface, 0.004 inch: time interval between dots, Ho second. 



R 

Figure 13.— Zone of amplification enclosed by neutral curve according to Schlichting. 
Theoretical total amplication from branch I to branch II as indicated. Experimental points 
from frequency of oscillations found on oscillograms. Branches I and II of neutral curve 
meet /?-axis at .R = ® . 

plification. The oscillations are therefore identified with 
those produced by instability of the boundary layer. 

2. THEORY OF BOUNDARY-LAYER OSCILLATIONS 

Evidence has been found of oscillations of rather high 
purity resulting mostly from random disturbances. The 
oscillations observed at a fixed point suggest the presence of 
a traveling wave in the boundary layer. A theoretical 
answer is sought for the question: How can random disturb- 
ances produce a wave that looks closely sinusoidal and has 
a single predominant frequency? 

The superposition of small disturbance velocities on a 
uniform velocity in one direction presents to a stationary 


observer the appearance of wave motion. The motion may 
be sinusoidal or irregular. For simplification, all the dis- 
turbance velocity components are assumed to be sinusoidal 
of the same frequency. In a uniform stream, this motion 
has the attributes of a traveling wave with a wave velocity 
equal to the mean velocity of the stream. The amplitude is 
the disturbance velocity with components u , v, and w , and 
the pressure involved is the dynamic pressure of the dis- 
turbance velocity. This kind of wave is involved in the 
application of stability theory to the boundary layer. Only 
two-dimensional motion has so far been successfully treated. 
In a nonuniform stream, such as a boundary layer, the wave 
velocity is not known at once but may reasonably be assumed 
(and can also be proved theoretically) to be less than the 
maximum velocity of the stream. The characteristics of the 
wave will depend on the characteristics of the stream, such 
as the velocity profile of a boundary layer, and on conditions 
imposed on the mean velocity and on the disturbance veloci- 
ties at the boundaries. The theory treats the influence of all 
these conditions and how they determine the wave length, 
the wave velocity, and the damping or the amplification of 
the wave. 

In general, the disturbances responsible for wave motion 
in a boundary layer are random in character and the result- 
ing wave motion is irregular. One of the important charac- 
teristics of the wave is its ability to maintain itself against 
the damping action of viscosity or even to grow by absorbing 
energy from the basic flow. Since this characteristic depends 
on the frequency, certain component frequencies of an irreg- 
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ular wave will be damped and others amplified as the wave 
travels downstream. It will be seen in section VI, 3, that 
this process accounts for the relatively high purity of the 
observed boundary-layer oscillations. 

Energy considerations alone form the basis of an elemen- 
tary treatment first used by Osborne Reynolds. According 
to the Reynolds concept, growth or decay of a disturbance 
depends on whether energy is transferred to the disturbance 
by absorption of energy from the basic flow or is extracted 
from the disturbance by the damping action of viscosity. 
The flow of energy then determines the stability of the 
boundary layer under the action of disturbances. Since the 
kind of disturbance must be specified, this method does not 
have general applicability. The flow of energy is, however, 
a fundamental physical concept involved in all theories even 
though not explicitly stated. As an introduction to the more 
advanced theories the energy relation will therefore be de- 
rived here. The development follows essentially that given 
by Prandtl in reference 20. 

(a) Basic equations and energy relation. — In general, 
only two-dimensional flow with two-dimensional disturbances 
has been treated successfully in stability theories. The basic 
flow is assumed to be steady and a function of y only. The 
basic flow in the boundary layer is then 

u=m 


v=o 



(4) 

(5) 


Equations (4) and (5) are the fundamental hydrodynamic 
equations for small disturbances on which are based all 
stability theories herein discussed. 

In order to derive an expression for energy balance, equa- 
tion (4) is multiplied by u and equation (5) by v and the 
resulting equations added to give the single disturbance 
equation 



( 6 ) 


It will be noted that the left-hand member of this equation 
gives the time rate of change of kinetic energy, due to the 
disturbances, of a particle moving with the basic flow. Each 
term on the right-hand side of equation (6) is now integrated 
over the region containing the disturbance, as explained in 
detail in reference 20. The resulting expression is 

w—'j J uv ^ ixdy -^ S(M~W dziy (7) 


The disturbances are a function of time as well as of x and y 
and are expressed by 

v=f2(x,y,t) 

The components of total velocity are therefore U+u and v. 
Since u and v are assumed to be small, the Navier-Stokes 
equations become, after squares and products of disturbance 
velocities and of their differential coefficients are neglected, 

du TT du dU ( di 2 U . d 2 u . d 2 u\ 1 / dP . dp\ 

c )t dx 1 dy v \dy 2 1 dx 2 ^~d?/ 2 / pydx^dx) 

( 1 ) 



where p is the pressure produced by the disturbances and P 
is the pressure due to the basic flow. 

The equation of continuity is 


By subtracting from equations (1) and (2), respectively, the 
corresponding Navier-Stokes equations for u=v=p= 0, 
the following two equations are obtained in terms of only 
the disturbances: 


where dE/dt now represents the time rate of change of kinetic 
energy, due to the disturbance, of the fluid within the region 
of integration. The first term on the right-hand side of 
equation (7) involves the shearing stress puv and represents 
the rate at which energy is absorbed from the basic flow. 
This energy is added to that already present in the dis- 
turbance. The second term on the right-hand side, which 
is obviously always negative, is the rate of dissipation of the 
energy of the disturbance by viscosity. The predominance 
of one effect over the other will determine the net change in 
the energy of the disturbance. 

The first term on the right-hand side of equation (7) has 
an interesting physical interpretation. First, if u and p 
differ in phase by 90°, as in the usual wave motion, uv will 
vanish over one complete cycle and no energy can be absorbed 
from the basic flow. Second, since dU/dy is positive, uv 
must be negative, if the disturbance is to receive energy 
from the basic flow. These and other conditions are treated 
in the more advanced theory to be considered next. 

(b) General disturbance equation and solutions. — Lord 
Rayleigh (reference 21) appears to have been the first to 
outline a more general mathematical theory to which numer- 
ous investigators have contributed. The mathematical 
difficulties are so great that the basic velocity profile ( U=f 
(y)) was usually approximated by a straight line or segments 
of straight lines. This imposed serious limitations on the 
generality of the results. Tietjens (reference 22) was the 
first to apply Rayleigh’s theory to profiles intended to repro- 
duce flow along a wall, but his solution suffered from the 
serious limitations just mentioned. Important advances 
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have been made by Tollmien (references 2 and 3) and 
Schlichting (references 4 to 7), who successfully applied the 
theory to curved profiles intended to represent actual velocity 
distributions. Their most complete solutions apply to an 
approximation of the Blasius distribution. For this dis- 
tribution Schlichting’s work (references 4 and 5) is par- 
ticularly complete and his results are used here for com- 
parison with experimental results. The following outline 
of the theory is taken largely from references 4 and 5 with 
occasional help from reference 2. The physical basis for the 
theory, the steps in the solution, and the more important 
results are given here. The mathematical aspects of the 
problem are discussed in the original references. Even in 
the original references details of the computations are lacking 
and checking of results is difficult. 

From the basic equations (4) and (5), the pressure terms 
are eliminated by differentiating equations (4) and (5) with 
respect to y and x, respectively, and subtracting the second 
from the first. The result is a linear homogeneous equation 
in u and v : 

& 2 u T d 2 u.Z>U()u b 2 U bvbU c) 2 v Ud 2 v _ 

dydt'~ ^ dxdy' dy dx ' v dy 2 ''dy dy dxdt dx 2 

v \dx 2 dy ‘ by 3 dx 3 dy 2 dx) 


The disturbance velocities may be expressed in terms of 
a stream function i/', such that 


u ~*y I 

V ~ dx J 


(9) 


(U—c) (F" — a 2 F) — U"F=^ ( F""-2a 2 F" + a i F ) (12) 

where the primes represent differentiation with respect to y. 
Equation (12) is a homogeneous linear differential equation 
of the fourth order with a general solution of the form 


F= C x Fi + C 2 F 2 + C 3 F 3 + C 4 F 4 (13) 

where F u F 2 , F 3 , and F 4 are the particular solutions and 
Cl, C 2 , C 8 , and C 4 are the constants of integration. The 
form of these particular solutions depends on the basic veloc- 
ity assumed, that is, on the boundary-layer profile. 

A solution of equation (12) was first obtained by Tollmien 
(reference 2) for a velocity distribution intended to approxi- 
mate closely a Blasius distribution. He considered the case 
in which amplification and damping were absent, that is, 
for real values of ft and c, and obtained the so-called “neutral” 
oscillations. Schlichting (reference 4) repeated Tollmien’s 
calculation and in addition treated the case for small ampli- 
fication. In a later paper (reference 5), Schlichting de- 
termined the distribution of amplitude of neutral oscillations 
across the boundary layer and investigated the energy 
balance. Schlichting^ work is therefore more complete 
and more suitable for comparison with experiment. 

Both Tollmien and Schlichting approximated a Blasius 
distribution by a straight line and a parabola as follows: 

For Og % ^0. 175, 

0 


For 0.175g | g 1.015, 


U 
V 0 


= 1.68 


Any periodic disturbance may be represented by a Fourier 
series. Since equation (8) is linear and homogeneous in 
u and v, its behavior may be investigated by using a single 
term of the series. This amounts to assuming a periodic 
disturbance whose stream function has the form: 

$=F(y) exp [i(ax~fit)] = F{y) exp [ia(x—ct)\ (10) 

where F(y) represents the initial amplitude of the stream 
function which depends only on y, a=2w/\ where X is the 
wave length, and t is the time. Since 0 and lienee c are 
generally complex quantities, equation (10) may be written 

$=F(y) exp {i[ax— (0 r +i0 4 )*]} 

=F{y) exp {ia[x— (c T +ic,)t]} (11) 

where 0 r , the real part of 0, is the angular velocity or 2wf 
(where / is frequency) ; 0„ the imaginary part of 0, is the 
coefficient of amplification or damping, depending on whether 
it is positive or negative; and c,=0 r /a is the velocity with 
which the phase progresses and will be termed simply 
“wave velocity.” 

When equation (8) is written in terms of the stream func- 
tion and substitution is made for in accordance with the 
second of equations (10), the result is 



For | >1.015, 



This gave a sufficiently close approximation to the Blasius 
distribution of U but was not a sufficiently good approxi- 
mation for U' and U". The latter quantities were therefore 
taken directly from the Blasius solution; only the following 
terms of the series, which obviously apply to small values of 
y/5, were used 



The displacement thickness 5*=J^ dy= 0.34:7 8. 

Schlichting gives 0.3415 for the integral and this value was 
used to compute the values of 5 in figure 21. 

It was assumed that U=f(y) and that F=0. This is 
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equivalent to assuming a nonthickening boundary layer, 
that is, constant <5, and is the basis for most of the adverse 
criticism of the theory. In a real boundary layer, 5 depends 
on x , TJ=j(x,y) r and WO; but the dependence of TJ on x is 
small compared with the dependence on y , and V is small 
compared with TJ. 

Since the purpose of the present paper is merely to in- 
dicate the method of solving the equation, the solutions will 
be indicated in general form. The calculations have been 
carried through to numerical results for only a few velocity 
profiles. 

Equation (12) is not readily solved in its entirety and 
certain approximations must be made that involve the fric- 
tion terms on the right-hand side of the equation. The 
simplest approximation is to neglect the friction terms en- 
tirely. In order to justify this, equation (12) is put in a 
dimensionless form that has the Reynolds number as one of 
the parameters. When the new independent variable y/5* 
is introduced and F=TJq5*<I>, the result is 

TV* <*"-2« ! **V'+«*s*V> 

(14) 

where the primes now denote differentiation with respect to 
y/5* and R is the Reynolds number. When R is sufficiently 
large, the right-hand side of equation (14) may be small 
enough to neglect. The effect of the friction terms is actu- 
ally negligible for all values of y/5 * except those near the 
surface and near a plane, 2 that is also near the surface, 
where the wave velocity is equal to the local stream velocity. 
It has been shown by Tollmien (reference 2) that c must 
always be less than U 0 ] a point will therefore always be found 
in the boundary layer where U=c. Exclusive of this point 
and the surface, two particular solutions 0! and 0 2 may be 
obtained by solving the frictionless equation 

(4>"-< x 2 8* 2 < t >) - 2=0 (15) 

The solutions are expressed as power series developed about 
y c , where y c is the position of the critical point, U=c. 
These solutions are 

♦.-V s * [■+«• ■ • ■] < 10 > 


(14), valid in the neighborhood of the critical point, results 
from neglecting all but the largest of the friction terms. 

For this purpose a new independent variable rj is intro- 
duced, defined by 

(18) 


From this it is seen that rj is a dimensionless variable which 
is equal to zero at the critical point y~y c , positive ioryJ>y c1 
and negative for ?/<Ct/ c Equation (14) is now written with 
rj as the independent variable for 0. Since this equation will 

U—c 

be used in the vicinity of y=y c , Tt~ can approximated 

U o 

/ jj ' \ 

i>y( er ) and TJ" by TJ v, c , where the primes again denote 
differentiation with respect to y/5*. The equation becomes 


i<f>" " (rj — 2ie 2 a 2 5 * 2 ) — 

4> (e jjr +>?eV5* 2 -ie 4 « 4 5* 4 )=0 (19) 

where 0 is differentiated with respect to tj. Since e is small, 
the terms containing e 2 and t 4 are omitted, and equation (19) 
becomes 

i<t>" " -\rrj<t>" =e jj-jr 0 (19a) 


Equation (19a) contains no singularity at U=c and a 
solution of this equation at the critical point indicates the 
form of 02 to be used when y—y c <i 0. This form is 


_ U" c ( 

0 2 — 'JJT <#>i 


log 


y-y c 

5* 




+ ^7T ) + 


7 i (y-y c ) , h (y-y c ) 2 

t — Ol 


(1 7a) 


By neglecting the term in equation (19a) containing e, 
the homogeneous equation 

i0"" + 770"=O (20) 


is obtained. The solution of equation (20) contains the 
two particular solutions 0 3 and 0 4 which are expressed in 
terms of Hankcl functions as follows: 


U" 

<t> 2 ^jpr<t>v log 


iy-y c ) A 

A* ’ T 


1 + ^1 


C y-yc ) 

5* 


b 2 


(y-y c ) 2 

5* 2 


(17) 


where U" c and U' c are the values of U" and TJ' at the 
critical point. While 0i is regular throughout the whole 
layer, 0 2 possesses a singularity at the critical point and a 
correction to equation (17) must be found by another ap- 
proximate solution that takes into account the friction at 
the critical point. This second approximation to equation 


<^ 3 . 4=2 dr > Jj dr, ( 21 ) 

where H (1),(2) are Hankel functions of the first and second 
kind. Equation (20) is valid all the way to the surface. It 
is found that 0 3 diminishes very rapidly from the surface 
outward and is therefore important only near the surface. 
On the other hand, 0 4 increases rapidly with y and will not 
satisfy the requirement that 0=0 at y= °° ; hence, (? 4 in 
equation (13) is set equal to zero. The general solution of 
equation (14) may then be written 


2 Since the theory applies only to a boundary layer of constant thickness, the problem is 
essentially one-dimensional and this plane is commonly referred to as the critical “point.’ f 


<t> = C\4> i + C 2 <#>2 + Cl 3 <j>3 


( 22 ) 
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(c) Characteristic-value problem. The stability investi- 
gation is a characteristic-value problem — that is, one in 
which the boundary conditions afford sufficient equations 
between the solutions at the boundaries to determine the 
values of the parameters for which equation (14) is satis- 
fied, those values being the characteristic values. At y=d 
where V= U 0 = Constant, the solutions of equation (14) are 
of the simple form 


Between these two equations any one of the three unknown 
parameters, ad*, It, and c/U 0 , may be eliminated and a 
relation between the other two found. This has been done 
by both Tollmien and Schlichting for c real, that is, for 
oscillations neither amplified nor damped. For this case, 
equation (27) may be written 

G n (rj 0n )=E n (a n d*,c r /U 0 ) (28) 


exp(-«5*|() 
exp («5* I*) 

The solution with the. positive exponent must be ignored as 
it is infinite at y=c°. As the outer boundary condition, 
then, 0'/ <t>= — ad*. Since 0 3 has already disappeared at the 
outer boundary, the result is 

C r i(0 , l 5 + «6*0 lg ) +C2(0 , 2 6 +«5*02 5 ) =0 
or 

C^ s +C 2 ^ 2d = 0 (23) 

where the sub-subscript 8 indicates the solution at the outer 
boundary and 

At the surface, y = 0 and 0=0' = 0. The surface boundary 
conditions are then simply 


where the subscript n denotes neutral oscillations. If c r is 
eliminated, a relation between ad * and R is found. This 
gives the theoretical neutral curve shown in figures 20, 25, 
and 27. By eliminating ad* the relation between c r /U 0 and R f 
shown in figures 26 and 33, is found. Recombination of 
these quantities results in P r vlU 0 2 , plotted in figures 13, 16, 17, 
19, and 24. 

The curves just mentioned extend only to the maximum 
values of R that are of practical interest. Tollmien points 
out in reference 2 that both branches I and II of these curves 
approach the It- axis as R approaches infinity. Branches I 
and II therefore meet at infinity and form a closed curve. 
Amplification occurs within the region enclosed by the curve. 
Damping occurs in all other regions. 

The values of the parameters are given in table III. The 
computations were made by a combination of analytical and 
numerical methods. It may be pointed out that the values 
obtained by Tollmien and Schlichting show some differences. 
These are probably due to differences in the details of the 
computations. Values plotted in the foregoing figures and 
appearing in table III are Schlichting’s results. 


Ci0i o + C 2 02 o + Cz<t> 3 0 = 0 (24) 

C Y 1 0'i o +r 2 0' 2o +C r 30 / 3 o=0 (25) 

where the sub-subscript 0 denotes solution on the surface. 

If 0 in equation (22) does not vanish identically, the 
following determinant of the boundary solutions, formed from 
equations (23), (24), and (25), must equal zero: 


*i. 0 

01 O <t> 2 0 <#> 3 0 

<t>'l 0 4>' 2 0 <t>'s 0 


= 0 


TABLE III 

THEORETICAL WAVE PARAMETERS FOR NEUTRAL OSCIL- 
LATIONS ACCORDING TO SCHLICHTING (REFERENCE 4) 




Branch I 


Branch II 

Cr 







Uo 

R 

ad* 

firV 

LV 

R 

a8* 

&rV 

Uo> 

0. 200 

7230 

0. 077 

2. 14X10-6 

37,600 

0. 149 

0. 79X10-« 

.250 

3010 

. 101 

8. 39 

12. 000 

. 188 

3. 92 

.300 

1530 

. 129 

25.3 

4,640 

.223 

14.4 

.325 

1150 

. 143 

40.4 

3.290 

.238 

23.5 

.350 

893 

. 159 

62. 3 

2,070 

.251 

42. 4 

. 375 

736 

. 181 

92.2 

1. 420 

.264 

69. 7 

.400 

633 

.205 

129.5 

1,020 

. 274 

107.4 

.420 

606 

.239 

165.7 

713 

.273 

161.0 


Solving the determinant gives 
fro 

frs 0 fr 0 frj — frofr^ 

The left-hand side of equation (26) is a function of rj 0 , 
where tj 0 is the value of rj at the surface. Since by equation 
(18) ry is a function of the parameters ad* and It, the left- 
hand side of equation (26) is also a function of these param- 
eters. The right-hand side of equation (26) is a function 
of the parameters ad* and c/U 0 . Denoting the left-hand 
side by G and the right-hand side by E gives 

GM=E(ad*,clUo) (27) 

Equation (27) is complex and may be separated into two 
real equations by separating real and imaginary parts. 


Schlichting extended the calculations to include small am- 
plification, applying only in the neighborhood of the neutral 
curve. The calculations were carried out, not by a direct 
solution of equation (27) but by a series development starting 
from the neutral equation (28). Thus G(rj 0 ) is found by 

G(rio) = G(rion) + (vo — Von)G / (rjon)+ . . . (29) 

and E(ad* , c/U 0 ) is found by 

E{ab*, C/U 0 ) =E n {a n b*, C T IU 0 ) + (a — a n ) 

(o-c r ) (^)+ • . • (30) 

With the values of G and E given by equations (29) and 
(30), amplifications were determined near and within the 
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space enclosed by the neutral curve ,#*= 0. Values beyond 
the range permitted by this method were determined by 
interpolation, by assuming a cubic equation for 0* at con- 
stant R that meets both branches of the neutral curve with 

the known slope ^ at these points. Schlich ting’s diagram 

showing the theoretical amplification is given in figures 27 
and 28. It will be noted that j 8 * has a maximum near the 
center of the zone and falls off toward both branches of the 
neutral curve. In accordance with equation ( 11 ), amplifica- 
tion of a disturbance is expressed by 

^?=exp £fi4t (31) 

where A x is the amplitude at t x and A 2 is the amplitude at t 2 . 
If t x is taken on branch I of the neutral curve and t 2 on 

branch II, the total amplification at constant ^ was calcu- 

Uo 

lated by Schlichting for several parts of the amplification 
zone. Sehlichting’s values are shown in figures 13, 19, and 
20 . The maximum total amplification up to a given value 
of r, when disturbances of all frequencies are present, is 
determined by the equation 

4M=exp (max J*™ P4t ) (32) 

Schlichting found that the relation between this maximum 
total amplification and the r-Reynolds number up to an R x 
of about 10 6 could be represented by the following interpola- 
tion formula: 


4^) = o.55 exp (0.555X10- 5 ^) (33) 

The solution thus far has determined several character- 
istics of the wave. For any given value of R and the fre- 
quency, it is known whether the wave will be amplified, 
damped, or neutral. For a given frequency the wave 
velocity is related to U 0 in a known way and is constant over 
the cross section of the layer. In like manner the wave 
length is known in relation to 6 *. 

(d) Distribution of amplitude, correlation, and energy 
balance. — The characteristic- value problem gives no infor- 
mation about the relative magnitudes of u and v nor how 
they are distributed through the boundary layer. How the 
boundary layer “vibrates,” whether as a whole or in parts, 
remains unknown. This is because the boundary conditions 
were used to determine only the parameters of the problem 
and not the constants of integration. 

In a second treatment of the theory applied to the flat 
plate (reference 5), Schlichting completed the investigation 
by determining the constants in equation ( 22 ) for neutral 
oscillations, that is, for j 8 *= 0 . He was thus able to calculate 
relative values of u and v and their distribution across the 
boundary layer. The physically important questions of the 
correlation between u and v and the energy balance were also 
answered. The same approximation to the Blasius distri- 
bution as used earlier was used here. 

8437 S4 — 50 3 


Since absolute values of u and v , which obviously depend 
on the intensity of the impressed disturbance, are of no 
interest in such a treatment, one of the constants was left 
undetermined and for convenience was made unity. Thus, 

1 

and, from the boundary conditions expressed by equations 
(23) to (25), the following values for C 2 and C z are obtained: 



When the values of the boundary solutions and their deriv- 
atives are known, C 2 and C z may be determined by equation 
(34). The general solution <f> is thereby determined except 
for an arbitrary intensity factor. From the nature of the 
particular solution <f> has real and imaginary parts <f> r and </>,-, 
which must be known with their derivatives in order to 
determine u and v. 

The stream function \p in equation (9) and the first of 


equations ( 10 ) are used to obtain the real parts of 77 and 

c'o 

v 

77 >since F=Z7 O 5*0, as follows: 

L'o 

u 

j T -=k \4>/ cos ( ax—f5 r t)—<j>i sin (ax—fi r t)] 

° [ (35) 

V 

jY=koi [<t> T sin (ax—/3 r t)A<l>i cos (ax—/3 r t)] 

Uo J 

where k is an arbitrary intensity factor. The corresponding 

root-mean-square values of u and v are 


Co _ -^V^ ,2+0/ 

V r ka / 


(36) 


(37) 


An expression for the kinetic energy is readily obtained by 
squaring equations (36) and (37) and adding. The correla- 
tion coefficient defined by 


becomes 




K = ~~r~i 

u'v 


<t>r' — 


V(«7 2 +<£/ 2 )(^ 2 +4n 2 ) 


(38) 


Values of u' / U 0 , v'/U 0 , and if were calculated by Schlichting 
for neutral oscillations corresponding to the following two 
points, one on each branch of the neutral curve: 


Branch 

R 

firV 

W 

ad* 

C r 

Uo 

I 

894 

62.3X10"« 

0. 159 

0. 35 

II 

2070 

42.4 

.251 

.35 


The value of k was chosen to give an average value of u'/U 0 = 
0.05 from 0 to 6. Table IV gives the values of u'/U 0 , v'/U 0l 
and K calculated by Schlichting. 

It will be seen in table IV that u' / U Q =0 at the surface and 
again at a point near y/8= 0.854 for branch I and near 
y/8= 0.774 for branch II, after which it again increases for 
greater y/8. The values of u'/U 0 are plotted in figure 21, 
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TABLE IV 

THEORETICAL NEUTRAL OSCILLATIONS TAKEN FROM 
SCHLICHTING (REFERENCE 5) 

[For average u'/U 0 =O.O5 from y/8 = 0 to y/8= 1] 


v 

, „ u' 

V 


b 

10 

t 0 

100 ~ 
Uo 

K 

Branch I; F , 894; ab *, *0.159 

0 

0 

0 

0 

.050 

.445 

.0546 

-.032 

. 090 

.647 

. 158 

-.061 

. 130 

.804 

.287 

-. 155 

. 170 

.922 

.434 

-. 169 

.209 

1.005 

.608 

-. 140 

.250 

1.016 

.795 

-. 109 

. 290 

1.003 

. 973 

-.064 

. 370 

.854 

1.322 

-.022 

.451 

.689 

1.605 

-.0075 

.531 

.524 

1.820 

0 

.012 

.371 

1.984 

0 

. 693 

. 244 

2. 098 

0 

.774 

. 121 

2. 170 

0 

.854 

.007 

2. 193 

0 

. 935 

.103 

2. 180 

0 

1.015 

. 211 

2. 108 

0 

1. 1 

.203 

2. 027 

0 

1.2 

. 193 

1. 935 

0 

1.3 

.185 

1.849 

0 

1. 4 

. 176 

1. 763 

0 

1.5 

. 168 

1.681 

0 

Branch II; F, 2070; <*5*, *0.251 

0 

0 

0 

0 

.029 

.661 

.0655 


.054 

.937 

.214 

-. 121 

.074 

1.038 

.359 

-. 170 

. 105 

1. 108 

.595 

-. 191 

. 157 

1. 127 

1.028 

-. 146 

.209 

1.086 

1. 445 

-. 076 

.250 

.998 

1. 755 

-.010 

. 290 

. 906 

2.04 

.015 

.370 

.686 

2. 50 

.019 

.451 

.504 

2. 87 

0 

..531 

.340 

3. 10 

0 

.612 

.205 

3. 26 

0 

. 693 

.087 

3.35 

0 

.774 

.022 

3. 37 

0 

.854 

.123 

3. 33 

0 

. 935 

. 216 

3. 22 

0 

1. 015 

.305 

3. 11 

0 

1. 1 

.292 

2. 925 

0 

1.2 

.272 

2. 72 

0 

1.3 

.252 

2. 53 

0 

1.4 

.234 

2.35 

0 

1.5 

.218 

2. 18 

0 


•$*=0.3415 


where the theory is compared with experiment. At the points 
within the boundary layer where u'/U Q vanishes, a 180° 
phase shift occurs in u , meaning that the instantaneous 
^-component of the oscillation velocity is decreasing on one 
side of this point while it is increasing on the other. No such 
phase shift occurs in v. Table IV also shows that correlation 
does exist between u and v and that uv is negative. This 
means that energy is drawn from the basic flow; and, in 
accordance with equation (7), this energy must be balanced 
by the dissipation if the oscillation is to be neither amplified 
nor damped. Sclilichting evaluated both terms in equation 
(7) by integrating from x=0 to x=\ and from y= 0 to y= co . 
The net flow of energy was found to be very nearly zero, 
which indicated that the oscillations were nearly neutral. 
This result was interpreted as proof of the correctness of 
the solution to the characteristic- value problem. 

A physical explanation of the correlation between u and 
v is attempted by Prandtl in reference 20. It appears that 
u and v would differ in phase by 90° and uv would vanish in 
consequence if there were no phase shift in v produced by the 
viscosity near the critical point. The viscosity is thus 
responsible for the flow of energy to, as well as from, the 
disturbance. 

(e) Solutions for other conditions. — Tollmien (reference 


3) has attempted to generalize the solution sufficiently to 
include boundary-layer profiles with an inflection point, 
such as occur with an adverse pressure gradient. For such 
profiles the sign of U" changes from positive to negative in 
passing outward from the surface and, at the inflection 
point, U" = 0. The velocity of the neutral wave is shown to 
be equal to the stream velocity at the inflection point, and 
an important consequence of this is the absence of a critical 
point at U=c r in the solution of the frictionless equation 
(15). The conclusion is that profiles with an inflection point 
are inherently unstable, although neutral oscillations can 
still exist. These considerations suggest that an adverse 
pressure gradient is favorable to amplification. 

Gortler (reference 23) has investigated the effect of wall 
curvature on the stability of the boundary layer and found 
that velocity profiles are unstable if a change of sign occurs 
in U" -\-U' /v, where r is the radius of curvature of the wall. 
If an inflection point (C7" = 0) exists, the point where 
U"+U'/r=0 is farther from the surface than the point 
where U"= 0 for a convex wall and nearer to the surface 
for a concave wall. If no inflection point exists, a point 
where ?7"+t77r=0 can occur only on a convex wall. 
It is concluded that convex curvature has a destabilizing 
effect and that concave curvature has a stabilizing effect. 

Solutions for profiles with inflection points are given in 
references 8 to 10. Rosenbrook (reference 8) investigated 
the flow through a diverging channel as the channel was 
towed through still water. Streamlines were made visible by 
fine particles of aluminum in the water and wavelike motions 
accompanied by the formation of eddies could be observed 
as well as the velocity distribution. A solution, following 
Tollmien in reference 3, was obtained for S-shape velocity 
profiles approximating those observed. Fair agreement was 
obtained between the observed and the computed wave 
lengths. Hollingdale (reference 9) obtained the solution 
appropriate to the laminar wake formed in one case behind a 
flat plate and in another behind an airfoil, both towed at low 
speeds through still water. Oscillating wakes were observed 
at certain Reynolds numbers and the wave length and the 
wave velocity were measured. With the appropriate inter- 
pretation of theoretical results, fair agreement between cal- 
culated and observed wave length and wave velocity was 
obtained. Savic (reference 10) obtained an exact solution of 
the frictionless equation (15) for the two-dimensional jet and 
obtained expressions for the neutral wave length and wave 
velocity in terms of characteristics of the jet. The theoretical 
results were compared with observations of wave length and 
wave velocity made on an acoustically sensitive jet and good 
agreement was obtained. 

Schlichting (reference 6) has determined the effect of 
centrifugal force by calculating the curves of neutral stability 
for a boundary layer developed on the inner wall of a rotating 
cylinder. In reference 7, Schlichting has determined the 
effect of density gradient on the stability of the boundary 
layer, using the approximation to the Blasius distribution 
given in part (b) of this section. Neutral curves outlining the 
zone of amplification are given. Stability is principally 
affected by a gravitational field in connection with the 
density gradient through the boundary layer. 
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The methods of Tollmien and Schlichting have been 
applied by Pretsch (reference 24) to the boundary layer on 
bodies of revolution. The approximate forms of the general 
disturbance equations from which the four particular solu- 
tions are obtained are found to be the same as those for the 
two-dimensional case; the solutions are therefore the same. 
The solutions are applied to the known velocity profile at the 
stagnation point of a sphere and the curve of neutral stability 
is computed. 

3. NATURALLY EXCITED OSCILLATIONS 

Naturally excited oscillations are those originating from 
the amplification of small disturbances naturally present in 
the stream. These oscillations were treated in section VI, 1, 
and their frequencies were experimentally obtained and 
plotted in figure 13. The concept of waves running down- 
stream through the boundary layer has been introduced; the 
behavior of the oscillations in relation to figure 13 can now 
be discussed. 

One important feature of figure 13 is that the values of 
amplification coefficient are positive within the zone enclosed 
by the neutral curve and have values reaching a maximum 
near the center and filling to zero on each boundary. Outside 
the zone p t is negative. All disturbances, whether damped 
or amplified, travel downstream and therefore in the direc- 
tion of increasing R since <5* increases with x. For example, 
a fixed speed U 0 and a disturbance of frequency such that 
P r v/U 0 2 in figure 13 equals 40X10 -6 are considered. Such 
a disturbance, if it began at the leading edge of the plate, 
would be damped as it passed along the boundary layer 
until it reached a point where II attained the value 1180. 
From this point to R= 2160 the disturbance would be am- 
plified and, as it passed beyond this point, would again be 
damped. If A x and A 2 are the amplitudes of the disturbance 
at the beginning and at the end of the amplification zone, 
the ratio A 2 tA x = 7±2. In other words, the amplitude lias 
increased 742 times while crossing the zone. This is the 
total amplification given by equation (31). A disturbance 
of any given frequency will then be most highly amplified 
when it reaches branch II of the neutral curve. 

Figure 8 shows that, in experiment, many frequencies are 
represented in the initial disturbance. The fact that the 
observed boundary-layer oscillations contain one predomi- 
nant frequency must be explained. In order to determine 
whether theory and experiment are in agreement, the place 
on the theoretical diagram where the predominant frequency 
should lie must be found. Instead of traversing the diagram 
of figure 13 with a wave of a single frequency as in the foregoing 
illustration, waves of many frequencies, covering a wide 
range of values of P^/Uq 2 , are supposed to enter the amplifi- 
cation zone through branch I and undergo amplification at a 
rate determined by their frequency as they proceed to the 
value of R at which the observation is made. Initial dis- 
turbances are, of course, present in the boundary layer at all 
values of R and disturbances enter the diagram at all points. 
Those entering this diagram ahead of branch I are reduced 
by damping below the prevailing level of the initial disturb- 
ances while those entering behind branch I cannot receive 
full amplification. Disturbances occurring at points corre- 


sponding to branch I have the largest initial value and also 
receive the highest amplification; therefore they are assumed 
to be the ones whose development should mainly determine 
the final result at the point of observation. It will be noted 
that amplification increases without limit as the frequency 
decreases if no limit is placed on the Reynolds number. 
The experimental observation, however, is made at a particu- 
lar value of R and hence of R x ] under this condition it is obvious 
that one particular value of B r vlU 0 2 will become more highly 
amplified than any of the others, this one having the maxi- 
mum total amplification given by equation (33). This fre- 
quency obviously is the predominant one that should be 
found on the oscillograph records. The remaining question 
now concerns the position of this most prominent frequency 
on the diagram, that is, where the experimental points should 
lie. The problem then becomes one of determining along 
what line, starting from branch I and running parallel to the 
7?-axis, the total amplification is a maximum up to a given 
value of R. By taking into account the values of Pi given by 
Schlichting it was found that, for uniform energy distribution 
with frequency, this line introduced into the diagram along 
branch I intersects the prescribed value of R somewhere 
between the center of the zone and branch II. For more 
accurate results, both the initial energy distribution and the 
damping prior to reaching boundary I must be considered. 
The experimental points should lie inside the zone near 
branch II. The agreement with theory is not perfect, since 
many points lie on branch II and outside the zone, but is as 
good as can be expected from observations of this sort. 

The process described here is believed to be one of con- 
siderable importance since it is the process by which laminar- 
boundary-layer oscillations develop from initially random 
disturbances. In order to develop a pure oscillation, initial 
disturbances must be so small that a large amount of am- 
plification can occur before transition destroys the laminar 
layer. 

It has been pointed out that the oscillations appear to be 
three-dimensional. This is to be expected because of the 
three-dimensional nature of the initial disturbances. The 
theory, however, concerns two-dimensional oscillations and 
how this should affect the agreement with experiment is not 
known. 

An apparent discrepancy between figures 7 and 13 must 
be pointed out. It will be observed in figure 13 that oscil- 
lations are indicated for R as high as 3120. This corre- 
sponds to an R x of 3.28 X10 6 and, according to figure 7, is 
in the transition region. The oscillograph records actually 
showed intermittent turbulence, but it was still possible to 
read oscillation frequencies from non turbulent parts of the 
records. The limits of the transition region are defined 
statistically by the surface-tube method. The discrepancy 
is therefore only apparent. 

In connection with the surface-tube method another ef- 
fect was noted — namely, that a rise in the surface-tube read- 
ing was noted a little ahead of the point where turbulence 
was indicated by the hot wire. This is believed to be 
caused by a slightly modified velocity distribution arising 
from the shearing stress given by uv in the oscillations. The 
lower boundary of the transition region in figure 7 is there- 
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fore a little too low, but the displacement is small and lias 
no significant effect on the shape of the curve. 

VII. ARTIFICIAL EXCITATION OF BOUND ARY-LAVEE 
OSCILLATIONS, ZERO PRESSURE GRADIENT 

It was soon realized that a study of boundary-layer oscil- 
lations could be carried out to better advantage if they were 
not caused by accidental disturbances occurring in the wind 
tunnel but were produced by a controlled disturbance of 
known amplitude and frequency at some chosen position. 
This was finally accomplished by the vibrating ribbon dis- 
cussed in section VII, 3. While a search was being made 
for devices of various kinds that would accomplish the de- 
sired result, effects of sound were tried. Interference from 
naturally excited oscillations was minimized by making the 
artificial disturbances much greater than the natural dis- 
turbances. The bolting-cloth damping screen was also in- 
stalled with the six wire screens, but the addition of this 
screen had no detectable effect on the natural oscillations. 

1. SOUND FROM LOUDSPEAKER 

The effect of sound intensity on the boundary layer is of 
interest because sound is an unavoidable source of disturb- 
ance both in wind tunnels and in free flight of power-driven 
aircraft. A 25-watt loudspeaker was accordingly installed 
in the top of the tunnel at the leading edge of the plate and 
fed by a variable-frequency oscillator through an amplifier. 
The intensity and the frequency were controllable. 

Interesting effects were easily demonstrated. For exam- 
ple, boundary-layer oscillations could be induced at will 
merely by choosing the right frequency for a particular 
position and speed. Transition could be moved 1 or 2 feet 
ahead of its natural position by the right combination of 
intensity and frequency. In general, a random noise from 
the loudspeaker produced similar results, but the effect on 
the oscillations was not so marked. These were casual ob- 
servations with no attempt at quantitative measurements. 

When quantitative work was attempted, it became ap- 
parent that the complicated sound field in the tunnel was 
a decided draw-back. Figures 14 and 15 show the distri- 
bution of sound intensity along the center of the working 
chamber 3 inches from the plate. Figure 14 shows a stand- 



leading edge of plate, measured along center 3 inches from plate. Uo= 0; frequency, 48 
cycles per second. 


ing wave at 48 cycles per second with either a node or a 
loop at the center of the plate, depending on the type of 
microphone — that is, crystal microphone, responding to 
sound pressure, or ribbon microphone, responding to particle 
velocity. This turned out to be the fundamental frequency 
of the 12-foot tube formed by the plate and the tunnel 
walls. The same standing wave could be set up without 
the loudspeaker when the exciting frequency came from the 
eight-blade wind-tunnel propeller running at a speed of 360 
rpm. Figure 15 shows the distribution at 60 cycles per 



Figure 15 .— Relative sound intensity from loudspeaker in ceiling of working chamber at 
leading edge of plate, measured along center 3 inches from plate. Frequency, 60 cycles per 
second. 


second. The ribbon microphone and the hot-wire anemom- 
eter show nearly the same distribution, as would be ex- 
pected. These are only a few examples of the sound field 
and the change in distribution with frequency that serve to 
complicate any studies with controlled sound intensity. 

When the loudspeaker was used to excite oscillations, two 
procedures were tried: (1) The sound frequencies that pro- 
duced oscillations of maximum amplitude at chosen velocities 
and positions were determined and (2) a sound frequency and 
position were chosen and the velocity was then varied until 
the oscillations were a maximum. By use of the sound fre- 
quencies read on the oscillator together with the measured 
U 0 and x , P T v/Uo 2 and R were calculated and plotted in figure 
16. The solid line is again Schlichting’s curve of neutral 
stability enclosing the zone of amplification. Since both 
procedures amounted to finding the condition for maximum 
total amplification, the points should lie within the zone near 



Figure 16. — Boundary-layer oscillations excited by loudspeaker. Points show sound fre 
quency giving maximum amplification. Curve outlines theoretical amplification zone. 
(See fig. 13.) 
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branch II to be in agreement with theory. In view of the 
difficulties encountered, the large scattering of the points is 
not surprising. 

Another application of the loud speaker consisted in di- 
rectly connecting the amplified output of the hot wire to the 
loudspeaker so that regeneration could occur. Disturbances 
picked up by the hot wire were thereby fed back as sound, 
which in turn reinforced the initial disturbance when the 
proper phase relation existed, causing the system to oscillate. 
The resulting frequency was read from a variable-frequency 
oscillator connected to the oscillograph to form Lissajous 
figures with the hot-wire signal. The frequency of regenera- 
tion was determined by the frequency response of the wire, 
the electrical circuit, the loudspeaker, and the tunnel as a 
whole, as well as by the preferred boundary-layer frequency. 
Alany frequencies could be found with no wind and only by 
rejecting these “false” frequencies could results of any relia- 
bility be obtained. The results are plotted on the usual 
diagram in figure 17. Alany of the points shown evidently 
are not related to the boundary-layer oscillations. 



R 

Figure 17.— Regeneration frequencies resulting from feedback from hot wire in boundary 
layer to loudspeaker. Curve outlines theoretical amplification zone. (See fig. 13.) 

2. SOUND THROUGH HOLE IN PLATE 

A localized disturbance was produced by bringing sound 
into the boundary layer through a hole in the plate. At a 
position 18 inches above the center line and 4 feet from the 
leading edge a %-incli hole was drilled through the plate and 
on the back was mounted a small headphone from a radio 
head set. This produced a three-dimensional disturbance 
emanating from a known location. 

A hot wire, mounted on the small sledlike head, was 
arranged for traversing upstream and downstream at a fixed 
distance from the surface, usually along a line slightly dis- 
placed from the hole. A survey of possible types of measure- 
ment showed that the method could be used only for obtain- 
ing the type of result shown in figure 18. With the re- 
generation hook-up, that is, when the amplified output of 
the wire fed into the headphone, the frequencies of oscillation 
were determined with the wire at various distances down- 
stream from the hole. The speed was constant during a run. 
Oscillation of the system occurred most easily for values of 
£7 0 less than 40 feet per second and with the hot wire less 
than 0.1 inch from the surface. The solid curves of figure 
18 show the discontinuous manner in which the frequency 
varied with the distance. Reversing the connections to the 



Figure 18. — Oscillations excited by sound through hole in plate. Points show regeneration 
frequencies resulting from feedback from hot wire in boundary layer to headphone. Hot 
wire moved fore and aft along line Me inch below r hole. Uo=23 feet per second. 


headphone changed the phase of the input by 180° and gave 
the results indicated by the broken curves. 

The explanation of this sort of performance is that proper 
phase relations had to exist for regeneration and that the 
frequency automatically decreased as the distance was in- 
creased in order to keep an integral number of waves between 
the hole and the wire. There existed a preferred frequency 
band, determined by the characteristics of the wire, electrical 
and acoustic circuits, and boundary layer. Jumps occurred 
when another wave between hole and wire was added to keep 
the frequency as near as possible to the center of the pre- 
ferred band. The distance between curves of one type along 
any line of constant frequency is thus the wave length 
corresponding to that frequency. 

It is believed that the frequency band was determined 
more by the boundary layer than by the rest of the circuit 
when the wire was far from the hole and that the preferred 
frequency was that for which amplification was a maximum 
over the path from hole to wire. The frequency midway 
between the extremes and as far as possible from the hole 
was therefore read from sets of curves for different wind 
velocities like the set in figure 18 and placed on the usual 
Schlichting diagram. The result is shown in figure 19. The 



Figure 19.— Middle frequencies of regeneration, wire far from hole, taken from curves of 
type shown in figure 18. Curve outlines theoretical amplification zone. Lines from points 
show length of amplification path. 
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path over which amplification could take place is indicated 
by the horizontal line connected to each point. 

Since the wave lengths corresponding to these frequencies 
were now known, it became possible for the first time to plot 
the results in terms of ad* and /?. The theoretical curve and 
the experimental results are shown in figure 20. The lines 
of total amplification for constant p r v/U 0 2 , corresponding to 
those in figure 19, are known. The paths over which ampli- 
fication could occur in the experiment are represented as be- 
fore by lines extending from each point. Figures 19 and 20 


<x6* 
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Figure 20.— Wave-length diagram showing zone of amplification according to Schlichting. 
Points show values of a5* for frequencies of figure 19. Lines from points show path of 
amplification. 

generally show the kind of results predicted by theory, if the 
frequency and wave length are determined mainly by the 
characteristics of the boundary layer. 

3. EXCITATION BY VIBRATING RIBBON 

(a) The ribbon. — All attempts to excite boundary-layer 
oscillations by wires vibrating in the boundary layer failed 
because the disturbance of the eddying wake from the vi- 
bration of wires large enough to affect the layer produced 
early transition. Flat strips outside the layer suffered from 
self-induced flutter and also had troublesome wakes that 
carried with them fluctuating pressure fields far downstream. 
A thin flat ribbon of phosphor bronze, placed edgewise to 
the flow in the boundary layer and only a few thousandths of 
an inch from the surface, was finally found capable of oscil- 
lating the boundary layer but produced no detectable effect 
when not vibrating. The ribbon was driven in and out from 
the surface by passing through it a small current of the de- 
sired frequency in the presence of a strong magnetic field 
from an electromagnet on the opposite side of the plate. 

Several ribbons were tried to find the best width, distance 
from surface, and length of vibrating span. All ribbons were 
0.002 inch thick. Widths ranged from Ke to ){ inch; 
distances from surface, from 0.006 to 0.015 inch; and vibrat- 
ing segments, from 6 to 12 inches. None of these factors 
appeared to affect the performance to any noticeable extent. 
A certain amount of tension was necessary for stability. 
When it was desired to place the resonant frequency above the 
working range, high tension was necessary. 

After this preliminary work, a ribbon 0.1 inch wide and 



3 feet long was installed 4 feet from the leading edge of the 
plate. At a distance of 6 inches from each side of the center 
line of the plate, two strips of Scotch cellulose tape were 
laid on the surface for insulation and spacing. Another strip 
of tape was then laid on the ribbon to hold it firmly in contact 
with the tape on the plate, except for the segment 12 inches 
long at the center. Rubber bands at the ends supplied the 
tension. The 12-inch segment at the center was 0.006 inch 
from the surface and was free to vibrate under the action of 
the current and the magnetic field. Examinations with a 
microscope showed that the ribbon vibrated to and from the 
surface in a single loop. The same oscillator and amplifier 
that had been used for the loudspeaker were found to furnish 
ample power for driving the ribbon. Different positions for 
the ribbon were chosen to cover a large range of Reynolds 
number but the method of attachment was always the same. 
Distance from the surface varied somewhat for different 
installations. 

(b) Procedure. — Two hot wires, one 0.010 inch and the 
other 0.110 inch from the surface, both sensitive only to u, 
were mounted on the small sledlike head and arranged for 
sliding upstream and downstream along the center line of 
the plate with the distance of the wires from the surface 
fixed. The fluctuations introduced into the boundary layer 
could be picked up by either hot wire. The wave form 
showed some distortion for an inch or so behind the ribbon 
but in general became nearly sinusoidal at greater distances. 
Frequencies amplified by the boundary layer showed increas- 
ing amplitude with increasing distance from the ribbon, 
while those damped showed decreasing amplitude. For 
much of the working region the oscillations could probably 
be regarded as two-dimensional but no attempt was made to 
verify this assumption. 

Just as for the loudspeaker and the headphone behind the 
hole in the plate, the system could be made to oscillate by 
connecting the amplified output of the hot wire to the ribbon. 
Curves like those shown in figure 18 were obtained and, from 
them, three points plotted in each of figures 19 and 20 were 
obtained. Since the ribbon could be used for applications 
that were believed to be more important, very little work was 
done with the regeneration hook-up. 

The wave length was determined in the following simple 
manner: The input to the ribbon from the oscillator was 
connected to one pair of plates in the cathode-ray oscillograph 
and the output from the wire was connected to the other pair. 
A stationary Lissajous figure consisting of a single closed 
loop was obtained since the frequencies of both input and 
output were the same. As the spacing between ribbon and 
hot wire was changed this figure changed from a straight line 
to an ellipse, then to a circle, again to an ellipse, and finally 
to a straight line inclined 90° to the first line. This indicated 
that the phase between input to the ribbon and output from 
the wire had changed by 180° and that the change in spacing 
between wire and ribbon was % wave length. With the 
oscillograph serving as the indicator, wave-length measure- 
ments were thus reduced to a measurement of the distance 
moved by the hot wire. The wave velocity was then ob- 
tained by multiplying the wave length by the frequency of 
the oscillations. 
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Determinations of amplification and damping were carried 
out by setting the ribbon into vibration at an arbitrary fixed 
amplitude and reading the mean-square hot-wire output 
voltage as the wire was moved away from the ribbon. The 
results obtained are given in section VII, 3 (d). The distri- 
bution across the boundary layer was determined by travers- 
ing normal to the surface at some fixed distance behind the 
ribbon with a single hot wire sensitive to u mounted on the 
traversing apparatus described in section III, 7. 

During the course of the work frequencies ranging from 10 
to 260 cycles per second were used. 

(c) Distribution of amplitude across boundary layer. — 
Since a theoretical distribution of amplitude across the 
boundary layer has been given by S chi ich ting (see sec. VI, 
2 (d), and table IV), the actual distribution is of interest as 
a test of agreement between theory and experiment. Tra- 
verses across the boundary layer were therefore made with a 
hot wire downstream from the ribbon with frequency, air 
velocity, and distance from the leading edge of the plate 
chosen to conform to the conditions for which the theoretical 
calculations were made. These conditions correspond to 
oscillations on the neutral curve specified as follows: 


Conditions required by theory 

Branch 

R 

w 


I 


894 



62.3X10-6 



II 


2070 



42.4 


Experimental conditions 




X 




Branch 

R 

firV 

w 

(ft) 


Uo 

(fps) 

Frequency 

(cps) 



Ribbon i 

Wire 




I 

902 

61.4X10-6 

0. 75 


1.09 

41.8 

103 

II 

2080 

40.6 

5.00 


5. 33 

46.0 

82 


i Ribbon about 0.009 in. from surface. 


The distributions are given in figure 21 with amplitude 
expressed in terms of u'/Uq . The solid curves show the 
theoretical distribution and the broken curves through the 
points show the observed distribution. In both cases the 
traverses were made 4 inches behind the ribbon. The vibra- 
tional amplitude of the ribbon was set arbitrarily and about 
the same for branches I and II. The smaller values of u'/U 0 
in branch II result from the smaller effect of the ribbon on 
the thicker boundary layer. Since the theoretical values 
given by Schlichting were based on an arbitrary average 
amplitude of u'/U 0 =0.05 from y = 0 to y=8 , Sclilichting’s 
values were reduced to give the average found by experi- 
ment; that is, areas under the theoretical curves were made 
equal to areas under the experimental curves from y = 0 to 
y=8. 

Although the general agreement between theory and ex- 
periment is good, perhaps the most striking confirmation of 
the theory is the 180° phase shift occurring at y/8= 0.7. On 
moving outward from the surface all oscillations disappeared 
at this point and farther out they reappeared but with phase 
opposite to those that had disappeared. This phenomenon 



Figure 21.— Distribution of amplitude of oscillations across boundary layer. Solid curves 
are theoretical according to Schlichting. Broken curves are experimental. 


is shown in a more striking manner by the oscillograms in 
figure 22. The pairs of traces are simultaneous records of 
oscillations picked up by two wires, one of which was at a 
fixed distance of 0.055 inch from the surface (lower trace) 
and the other at varying distances, beginning at 0.080 inch 
(upper trace). The simultaneous records were made pos- 
sible by an electronic switch that impressed the separately 
amplified voltages of the two wires on the oscillograph, 
alternately and in rapid succession. The switching from 
one to the other accounts for the dotted appearance of the 
trace. The waves are in phase for the first four records; the 
amplitude of the upper wave decreases with each succeeding 
record to zero in the fifth. Beginning with the sixth, an 
oscillation reappears in the upper trace but 180° out of phase 
with that in the lower trace. The oscillations in these 
records lie near the center of the amplification zone at 
72=1720. 

The type of distribution shown in figure 21, including the 
phase reversal, was found for amplified, damped, and neutral 
frequencies and for oscillations as close as 1 inch behind the 
ribbon and as far downstream as the oscillations remained 
free from transition disturbances. The fact that the main 
features of the distribution were found 1 inch behind the 
ribbon shows that the forced oscillations in the boundary 
layer very soon reach a steady state. 
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Figure 23.— Growth and decay of u-component of oscillations produced by vibrating ribbon 
4 feet from leading edge of plate. To, 2 inches downstream from ribbon; Uq = 64 feet per 
second. 


denotes a position 2 inches downstream from the ribbon. 
The ratio u'/u 0 ' is the square root of the meter reading at x 
divided by the square root of the meter reading at x 0 . 
According to equation (31) 


or 


u' ffO) 

57 ““ p J.w 


&idt 



Pidt 


Differentiating with respect to t and using the relation 
dx/dt=c T , where c, is the wave velocity, give 


Pi= 2.3 



(39) 



Figure 22.— Oscillograms showing phase reversal in u-component of oscillations excited by 
vibrating ribbon. Simultaneous records from two hot wires located 1 foot downstream 
from ribbon. Lower trace, hot wire 0.055 inch from surface; upper trace, hot wire at various 
distances from surface; ribbon 3 feet from leading edge of plate; frequency, 70 cycles per 
second; Uo= 42 feet per second. 

(d) Amplification and damping/ — When the wire was 
kept at a fixed distance from the surface and moved down- 
stream from the ribbon, the amplitude of the forced oscilla- 
tion either increased, decreased, or remained the same, de- 
pending on the frequency, the velocity, and the Reynolds 
number. The amount of the change with distance could be 
determined by reading the mean-square value of the fluctu- 
ating voltage across the wire on the output meter of the 
turbulence amplifier. In order to eliminate as far as pos- 
sible the extraneous disturbances introduced by turbulence 
and by the noise of the tunnel, the signal from the wire was 
passed through an electrical filter, which w T as tuned in each 
case to pass only a 5-cycle band around the frequency being 
measured. The results of a typical run with the ribbon at 
the 4-foot position are shown in figure 23. In this figure £o 


The w^ave velocity is obtained by multiplying the measured 
wave length by the frequency. If the changes in u' lu 0 ' with 
distance shown in figure 23 are truly effects of amplification 
and damping, the slopes of the curves in this figure serve to 
determine 0*. 

Before these slopes can be used with confidence, effects 
introduced by the thickening of the layer with increasing x 
and effects from possible changes in the distribution across 
the layer also with x must be investigated. Errors can arise 
from these sources because of the experimental procedure — 
namely, the determination of growth or decay of amplitude 
by varying the distance between the ribbon and the wire 
while keeping the wire at a fixed distance from the surface. 
From the evidence in section VII, 3 (c), changes in the dis- 
tribution are small and can cause but little error; however, 
the increase of 5 with x may give rise to an apparent amplifi- 
cation or damping, depending on the position of the wire 
in the layer. In all cases in which data on amplification 
and damping were obtained, the wire was in the rising part 
of the distribution, that is, values of y/b below those for 
maximum u'/Uq in figure 21. Consequently, the increasing 
8 with increasing x introduced an apparent damping, which 
caused all observed amplification to be less than the true 
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value and all observed damping to be too great. If it is 
assumed that the inner part of the distribution curve is a 
straight line passing through the origin, the magnitude of 
the error in 0* calculated from equation (39) by use of ob- 

\ Cf 

served slopes is —• Any curvature in the distribution 

Z Ju 

curve like that shown between zero and the peak value of 

1 c 

u' / U 0 in figure 21 lessens the error. Since ^ ^ therefore 

represents an upper limit and since its value turned out to 
be about the order of the random deviation in results, fit was 
calculated by equation (39) with the use of the experimentally 
determined slopes without correction. 

In figure 23 the curve for 40 cycles per second shows damp- 
ing tlu’oughout. As the frequency is increased amplification 
takes place, reaches a maximum at 120 cycles per second, 
and then decreases until 180 cycles per second again shows 
damping. In the range of Reynolds number represented 
by this range in x, all frequencies below about 50 and above 
180 cycles per second are damped, while a band of ampli- 
fied frequencies lies in between. 

Of considerable interest are the frequencies for which 
there is neither amplification nor damping and the corre- 
sponding values of R. This amounts to finding points of 
zero slope on curves like those shown in figure 23. Such 
points were obtained from a number of runs at various 
wind speeds with the ribbon at the 4-foot position and 
again with the ribbon 9 inches from the leading edge. Values 
of p r v/U 0 2 and the corresponding values of R were computed 
and plotted in figure 24. The solid curve is again Schlicht- 
ing’s neutral curve and the broken curves are those defined by 
the experimentally determined points. The agreement is 
very satisfactory, except for points above £br/?7o 2 =200X 10“ 6 . 
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Figure 24. — Frequency of neutral oscillations excited in boundary layer by vibrating ribbon. 
Solid curve is neutral curve according to Schlichting. Broken curves are neutral curves 
defined by experimental points. 


By use of the measured wave lengths, a8* was calculated 
for the neutral frequencies and the results were plotted on 
the wave-length diagram shown in figure 25. Schlichting’s 



R 


Figure 25.— Wave length of neutral oscillations excited in boundary layer by vibrating 
ribbon. Theoretical neutral curves shown solid, neutral curves defined by experimental 
points shown broken. Closed circles, branch I. Open circles, branch II. 

neutral curve is shown by a solid line and the broken 
curves are those defined by the experimental points. 3 * * Open 
circles denote points on branch II, and closed circles, points 
on branch I. 

Still another type of representation involving both fre- 
quency and wave length is the diagram c r /Uo and R in figure 
26. The points and the broken curves show the wave veloci- 
ties found for the neutral oscillations. The solid curve is 
Schlichting’s theoretical curve. 

It is again pointed out that both branches of the theo- 
retical curves in figures 24 to 26 meet the R - axis at R= °° 



Figure 26.— Wave velocity of neutral oscillations excited in boundary layer by vibrating 
ribbon. Solid curve is neutral curve according to Schlichting. Broken curves are neutral 
curves defined by experimental points. Closed circles, branch I. Open circles, branch II. 

and thus extend beyond the region indicated. Transition 
prevented experimental determinations above 7? = 2800. 

The wave-length diagram is repeated in figure 27 and this 
time shows, in addition to the theoretical neutral curve, the 

3 The theoretical neutral curve given by C. C. Lin has been added in figure 25. Lin’s 

curve was not available when this work was originally written. The curve is taken from C. 

C. Lin, “On the Stability of Two-Dimensional Parallel Flows," Quart. Appl. Math., vol. 

Ill, no. 2, July 1945, pp. 1 17-142; no. 3, Oct. 1945, pp. 218-234; and no. 4, Jan. 1946, pp. 277-301. 
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Figure 27. — Contours of equal amplification according to Schlichting. Values of ^,-5*/Ub 
(all values to be multiplied by 10 -4 ) opposite points are amplifications determined by 
experiment. Faired experimental contour of zero amplification shown by broken curves. 


curves of equal amplification within the zone. These con- 
tours were copied directly from a figure given by Schlichting 
in reference 4. The broken curves are defined by the neutral 
points, which are not shown in this figure. The points with 
the values opposite (all to be multiplied by 10 -4 ) are the 
observed values of Pi8*IU 0 obtained from values of the slopes 
other than zero. Negative values signify damping. 

The distribution across the amplification zone at particular 
values of R is shown in figure 28. The theoretical curves, 
drawn solid, are plotted from values read from Schlich ting’s 
figure in reference 4. The broken curves are drawn through 
the experimental points. 



Figure 28.— Distribution of amplification across zones of figure 27. Solid curves are 
theoretical. Broken curves are experimental. 


The comparison between theory and experiment given in 
figures 24 to 28 shows, on the whole, better agreement than 
might be expected when it is considered that at the outset 
the very existence of such phenomena had not been proved. 
The observed amplification shown in figures 27 and 28 is gen- 
erally less than the theoretical, but it will be recalled that the 
theoretical calculation yielded accurate values only near the 
neutral curve. In the neighborhood of the low Reynolds 
number limit the experimental inaccuracies are large because 
of the small amplification and damping. The maximum 
values of /3 r r/C7 0 2 , &8*, and c r /U 0 and the minimum value of 
R for which amplification can occur as determined approxi- 
mately by experiment are compared with the theoretical 
limits in the following table: 



Experiment 

(approx.) 

Theory 

Schlichting 

Tollmien 

R 

450 

575 

420 

0tV 

Uo* 

400X10- 6 

178X10-6 


a8* 

0.4 

0. 278 

0. 367 

Cr 

Uo 

0. 43 

About 0.42 

0. 425 


Tollmien’s limiting values, on the whole, agree better with 
experiment than Schlichting^ values. A comparison of the 
experimental results in figure 25 with Tollmien's diagram of 
a<5* plotted against logi 0 R in reference 2 shows better 
agreement in general than with the corresponding diagram 
given by Schlichting in figure 25. In spite of this, Schlicht- 
ing’s results were used because they were more complete 
than Tollmien’s. 

VIII. EFFECT OF PRESSURE GRADIENT ON BOUNDARY- 
LAYER OSCILLATIONS 

Figure 29, obtained in connection with the earlier work 
presented in section VI, 1, shows the effect of pressure gra- 
dient on naturally excited oscillations. The pressure dis- 
tribution along a part of the plate is shown to the left in this 
figure. The reduction of amplitude in the region of falling 
pressure is very marked. Accelerated growth of amplitude 
is evident in the region of rising pressure. It was accordingly 
inferred that negative pressure gradients decreased amplifi- 
cation (or increased damping) of the oscillations while posi- 
tive gradients increased amplification. A quantitative in- 
vestigation of this effect was therefore undertaken with 
oscillations forced on the boundary layer by the vibrating 
ribbon. 

When the ribbon was in the 3-foot position, the regular 
procedure of measuring amplifications and damping and of 
mapping the neutral curves was carried out for pressure 
distributions B, C, D, and E shown in figure 30. The cor- 
responding pressure gradients are given in table V. 

In figure 30, q x is the dynamic pressure at any distance x 
from the leading edge while q 0 is a convenient reference 
pressure. Curve A shows one of the distributions for the 
condition generally termed “zero pressure gradient.” This 
condition was obtained by giving sufficient divergence to the 
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Figure 29.— Oscillograms of w-fluctuations showing effect of pressure gradient on boundary- 
layer oscillations. Distance from surface, 0.021 inch; Uo= 95 feet per second; time interval 
bet ween dots, Ho second. 


TABLE V 

PRESSURE GRADIENTS 

\?i Ax/ 

[qi, dynamic pressure outside boundary layer at 3.33 ft from leading 
edge of plate; P , pressure near surface at distance x from leading 
edge] 


Distribution 

Average gradient 
from 1 to 6 feet 
(percent) 

Local gradient 
(percent) 

A 

-0. 11 

—0.78 from 4 to 5 ft 

B 

1.62 

2.48 from 3 to 3.67 ft 

C 

-2. 38 

-1.36 from 3 to 3.67 ft 

D 


12.0 from 3 to 3.67 ft 

E 


—9.8 from 3 to 3.67 ft 


1 Note that pressure gradient and dynamic pressure gradient have opposite signs. (See 
fig. 30.) 

adjustable side walls of the working chamber to overcome the 
normal pressure drop. The shape of the curve in the first 6 
inches was determined by the flow around the leading edge of 
the plate and hence could be changed but little by the setting 
of the walls. The walls also had little effect on the small 
local variations, since these were found to be produced by 
slight waviness of the surface. Curve A is thus a condition of 
zero gradient only when the length of the interval and the end 
points are properly chosen. With a pressure distribution 
given by curve A, the experimental values of velocity dis- 
tribution and 6* were in good agreement with the Blasius 
values; hence, the approximation to zero gradient is probably 
sufficiently good. 

Curve B, which shows the pressure rising along the plate, 
was obtained by giving the walls their greatest possible 
divergence. Curve C, which shows the pressure falling, was 



obtained with the walls forming a converging passage. The 
steep positive gradient shown by curve D was produced by a 
vertical airfoil placed 10 inches from the side of the plate 
running from the lloor to the ceiling of the tunnel, and the 
steep negative gradient of curve E was produced by blocking 
the stream with a plywood panel % inch thick and 28 inches 
wide extending from the floor to the ceiling placed parallel to 
and 10 inches from the plate. 

The results obtained on the effect of pressure gradient are 
given in figures 31 to 34 for distributions B and C and will 
be discussed in detail later in this section. Corresponding 
figures for distributions D and E could not be obtained 
because of the large effect of these gradients on the oscilla- 
tions. With the large gradients, the general impression 
given by observing the oscillations was that amplification 
always occurred in gradient D regardless of the frequency 
and Reynolds number, while damping always occurred in 
gradient E; however, this statement must be somewhat 
qualified. For example, in E fluctuations of all frequencies 
were damped for values of R up to 2600, which was the 
highest value obtainable with the ribbon at the 3-foot posi- 
tion. It is quite possible that a zone of amplification lay 
beyond the range of the experiment. In the case of gradient 
D, the investigation had to be restricted to low Reynolds 
numbers because of the occurrence of transition. The few 
observations made here indicated a possible low limit of 


26 


REPORT NO. 90 9 NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 


Reynolds number below which damping occurred, but the 
results are in doubt because of a possible scale effect on the 
airfoil used to produce the pressure gradient. The pressure 
rise may have been lessened by laminar separation on the 
airfoil at the low velocities and the positive gradient may 
have become less than that given by D. Accurate measure- 
ment of pressure distribution at the low velocities was 
difficult. 

Before presenting the results for the two small gradients 
B and C, attention must be called to the effect of pressure 
gradient on <5*. Since <5* enters into certain dimensionless 
quantities used in plotting the results, its actual value must 
be known for each case. The value of 8 * was measured only 
for distribution B and the result was 


Observed 5* 
Blasius 8 * 


= 1.09 at 2=3.33 feet 


For distribution C, 8 * was estimated by the aid of the 
K&rmdn-Pohlhausen method by assuming a linear pressure 
distribution. The validity of the method was tested by 
making a similar calculation for B and obtaining good agree- 
ment with the measured 8 *. The estimated result for C was 


Calculated 8 * 
Blasius 8 * 


= 0.93 at £=3.33 feet 


The values of 8 * for the Blasius distribution were multiplied 
by these factors to give the value of 6* pertaining to gradients 
B and C. 



Figure 31.— Frequency diagram of neutral oscillations. Experimental curve for zero pres- 
sure gradient shown without points. Positive gradient from distribution B and negative 
gradient from distribution C of table V. 


Figure 31 is the much-used frequency diagram showing 
the results for zero, positive, and negative gradients of dis- 
tributions A, B, and C, respectively. The solid curve is the 
same as the broken curve of figure 24 for zero gradient, 
shown here without the experimental points. The points for 
negative gradient lie close to this curve; the points for posi- 
tive gradient generally lie outside the zone of amplification 
for zero gradient. The positive gradient has therefore 
widened the amplification zone as expected. Contrary to 
expectation, the negative gradient failed to reduce the size 
of the zone. No observations were made below 72=800 . 

Since wave lengths of the oscillations were measured, the 
data for these same frequencies are given in a wave-length 
diagram (fig. 32) and a wave-velocity diagram (fig. 33). It 



R 


Figure 32.— Wave-length diagram of neutral oscillations. Experimental curve for zero pres- 
sure gradient shown without points. Positive gradient from distribution B and negative 
gradient from distribution C of table V. 



R 


Figure 33.— Wave-velocity diagram of neutral oscillations. Experimental curve for zero 
pressure gradient shown without points. Positive gradient from distribution B and 
negative gradient from distribution C of table V. 
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will be noted that figure 32 is unlike the usual wave-length 
diagram (for example, fig. 25) in that the ordinate is now 
av/Uo instead of ad *. This change was made in order to 
eliminate 5* from one of the coordinates and thereby to re- 
strict errors occurring in the estimated value of 8* with pres- 
sure gradient to the Reynolds number coordinate. The new 
ordinate is simply c*5* divided by the Reynolds number, 
that is, 

o:5* av 


ber (see fig. 27), the difference between 77=1840 for zero 
gradient and 77=1810 for negative gradient may account 
for the reduced amplification with negative gradient. The 
failure to show conclusively an effect of the small negative 
gradient may be due in part to the small local value of the 
gradient shown in table V. The difference between 77 =2200 
for zero gradient and 77=2050 for positive gradient is sig- 
nificant, since 7? =2050 for zero gradient would lie below 
7?=2200. 

IX. TRANSITION FROM OSCILLATIONS TO TURBULENCE, 
ZERO PRESSURE GRADIENT 


In both figures 32 and 33, the solid curves are the experi- 
mentally determined neutral curves for zero gradient with 
the experimental points omitted. 

Figures 32 and 33, like figure 31, show an effect of the 

positive gradient but little or no effect of the negative 

gradient. These figure show only the neutral points and 

so merely define the boundaries of the amplification zone. 

The effect of pressure gradient on the amplification within 

the zone is demonstrated in figure 34. This figure, which is 

similar to figure 28 except that zones are from the frequency 

diagram instead of the wave-length diagram, shows the 

( 3 d*\ 
actually jj~j 

across the zones of figures 31 at certain fixed values of the 
Reynolds number. The positive and negative pressure 
gradients signify as before pressure distributions B and C, 
respectively. For comparison, the experimental curves for 
zero gradient are given as solid lines. The comparison is 
best made between pairs of curves having Reynolds numbers 

Bid* 

most nearly alike. Since yj- increases with Reynolds num- 



Figure 34.— Distribut ion of amplification across zones of figure 31, showing effect of pressure 
gradient. Positive gradient from pressure distribution B and negative gradient from dis- 
tribution C of table V. 


The distinguishing feature of transition with low stream 
turbulence is the occurrence of boundary-layer oscillations 
that grow until the layer becomes turbulent. Under such 
conditions, transition to turbulent flow is simply transition 
from a particular type of oscillatory motion to turbulent 
motion. The oscillations are now understood and the final 
step in an understanding of the origin of turbulence involves 
an investigation of this transition process. 

Although transition in highly disturbed flow does not fall 
into this category, the two are similar in that both have to 
do with the effect of a disturbance large enough to bring 
about turbulent flow. The effect, however, depends on the 
type of the disturbance as well as its magnitude; and, as 
has been shown, an amplified disturbance does not have the 
random character of an initially large disturbance. Transi- 
tion from an amplified disturbance promises to be of con- 
siderable importance in modern low-turbulence wind tunnels 
and in the atmosphere. 

The vibrating ribbon to produce oscillations and the hot- 
wire apparatus to examine their development were suitable 
tools for the investigation of transition. The ribbon was- 
placed at the 3-foot position and two hot wires at different 
distances from the surface were arranged for traversing par- 
allel to the plate through the transition region. The Reynolds 
number was always chosen so that the undisturbed boundary 
layer (ribbon not vibrating) was laminar at the hot wires. 
Representative samples of the results obtained are shown in 
figures 35 to 37. 

Figure 35 shows what happens at a fixed point 2.5 feet 
downstream from the ribbon as the vibrational amplitude of 
the ribbon is varied. The pairs of traces show simultane- 
ously the response of the two wires. The values of gain 
control setting on driving amplifier are roughly proportional 
to ribbon amplitude. The amplitude of the waves on the 
traces cannot be quantitatively related to fluctuations of 
velocity because the sensitivity of the wires was unknown. 
The sensitivities of the two wires were different, and so com- 
parison of amplitudes at inner and outer wires was impos- 
sible. It is known, however, that velocity increases produce 
upward displacement of the traces and that velocity de- 
creases produce downward displacement. The direction of 
time in this figure is from left to right. 

This figure shows that the boundary layer at a given point 
can be changed from laminar to turbulent by varying ribbon 
amplitude and shows qualitatively the progressive nature of 
the change. Examination of successive pairs of traces indi- 
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Figure 35.— Oscillograms of ^-fluctuations showing development of turbulence from oscilla- 
tions produced by vibrating ribbon. Simultaneous records from two hot wires located 
2.5 feet downstream from ribbon. Upper trace, hot wire 0.014 inch from surface; lower 
trace, hot wire 0.114 inch from surface; ribbon 3 feet from leading edge of plate; frequency, 
70 cycles per second; Uo=50 feet per second. 


cates that fairly regular waves first make their appearance, 
being first in phase at the two wires and later in the fifth 
and sixth records, in and out of phase in random fashion as 
the amplitude increases. The distances from the surface are 
such that an in-phase condition would be expected from the 
distribution in figure 21. A growing distortion of the waves 
then appears until, in the seventh record, bursts of high- 
frequency fluctuations occur. The latter, identified by the 
jagged appearance of the trace, are bursts of turbulence. 
The final record shows the completely turbulent state. The 
dotted character of the traces, due to the electronic switch, 
obliterates details of rapid fluctuations. A similar succession 
of changes is shown in figures 36 and 37. 

In figure 36 simultaneous records are again shown, this 
time from two wires on the same mounting moved down- 
stream from the ribbon while the amplitude of the ribbon 
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Figure 36.— Oscillograms of u-fluctuations showing development of turbulence from oscilla- 
tions produced by vibrating ribbon. Simultaneous records from two hot wires. Upper 
trace, hot wire 0.011 inch from surface; lower trace, hot wire 0.112 inch from surface; ribbon 
3 feet from leading edge of plate; frequency, 70 cycles per second; Uo=50 feet per second. 


remained fixed. Figure 37 shows the same condition but 
for the two traces photographed separately in order to 
eliminate the dotting effect of the electronic switch. The 
difference in the appearance of the records in these two 
figures is due more to the different magnification than to 
the character of the traces. Since the frequency of 70 
cycles per second is amplified by the boundary layer, the 
amplitude increases downstream and increasing distance is 
in one respect like increasing ribbon amplitude in figure 35. 
The two cases are not equivalent, however, because the 
boundary-layer thickness increases from record to record 
in figures 36 and 37. 

From calibration data, obtained in connection with 
figures 36 and 37, it was possible to relate two particular 
displacements on the record to a percentage change in the 
mean velocity at the position of the wire. These two per- 
centages are indicated at the right of each trace, an upward 
displacement indicating an increase and a downward dis- 
placement a decrease from the mean local velocity at the 
hot wire. The percentages are the same for all records in 
figure 36. 

For fluctuations as large as most of those shown, it is not 
possible to indicate the zero- or mean-velocity position of 
the trace. This is because of the nonlinear characteristics 
of the hot wire which produces greater displacements of the 
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Figure 37.— Oscillograms of ^-fluctuations showing development of turbulence from oscillations produced by vibrating ribbon placed 3 feet from leading edge of plate. Freq uency, 70 

cycles per second; (7o=50 feet per second. 


trace for velocity decreases than for increases. A certain 
amount of distortion and some error in the records will 
result from this effect. It is certain, however, that crests 
are above the mean velocity while troughs are below the 
mean and that the distance between the two is approximately 
the double amplitude. The scale at the right furnishes a 
measure of the double amplitude. 

In the succession of changes occurring from point to 
point through the transition region, the oscillations first 
increase in amplitude and, once large, undergo increasing 
amounts of distortion. Turbulence appears at first in short 
intervals of the trace and finally for the full length of the 
trace. Randomness occurs gradually, and this suggests a 
continuous transition from wave to turbulent motion. It is 
not difficult to imagine a process here like that often assumed 
for the formation of eddies from a free vortex sheet. The 
sheet is imagined to take first a wavelike character, then as 
the wave grows to curl up into discrete eddies. The dis- 
turbed laminar boundary layer may be regarded as a wavy 
vortex layer with the wave progressively increasing in ampli- 
tude and distorting until discrete eddies are formed. The 
eddies themselves are unstable and soon break up into a 
diffusive type of motion which characterizes turbulent flow. 
If the validity of this picture is granted, transition is under- 
stood when the wave motion is understood. An under- 
standing of transition is therefore close but only close because 
the wave motion is understood only when the amplitude is 
small. Furthermore, the important question concerning 
where along the surface eddies first form is still unanswered. 


Another process may also be imagined. It will be noted in 
figures 35 to 37 that transition takes place in a more orderly 
manner at the greater distance from the surface. Here 
turbulence appears in certain parts of the wave cycle, gen- 
erally in the sharp-pointed low-velocity part of the cycle. 
This behavior lends support to the view held by some investi- 
gators that transition occurs when the flow comes to rest or 
reverses direction near the surface (separation point). While 
the amplitudes shown do not reach the necessary 100 percent, 
it seems possible that this value would be found at some point 
nearer the surface. Against this argument is the fact that 
the records made nearer the surface do not show the turbulent 
bursts confined to the low-velocity side of the cycle. The 
argument is further weakened by the possibility of accounting 
for the observed effects at the greater ^/-values by the differ- 
ence between the distribution of mean velocity characteristic 
of a laminar layer and that characteristic of a turbulent layer. 
This is to say, that the observed effect is caused by the turbu- 
lence rather than a cause of the turbulence. 

The idea of transition by separation has been formulated 
by Taylor (reference 1) by assuming intermittent separation 
or near separation to be caused by the local pressure gradients 
accompanying disturbances. The assumption is that separa- 
tion is imminent when the Karman-Pohlhausen parameter A, 
derived from pressure gradients either steady or fluctuating, 
reaches a certain negative value. A steady adverse pressure 
gradient large enough to cause laminar separation is known 
to bring about transition rather than complete separation 
when the Reynolds number is sufficiently high. By analogy 
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intermittent separation is supposed to have the same effect. 
For isotropic turbulence in the surrounding flow the fluctuat- 
ing pressure gradient is known in terms of measurable quan- 
tities, namely, u' and the scale. Taylor’s relation is 



where L is the scale and F{R X ) expresses a function of the 
^-Reynolds number of transition. 

In essence Taylor’s relation rests on the concept of velocity 
fluctuations impressed on the boundary layer from pressures 
outside. When velocity fluctuations belong to the natural 
wave in the boundary layer, Taylor’s relation can have 
little bearing on the problem. There is therefore no theory 
connecting oscillations to transition. The idea of transition 
by separation leads to the concept that transition results 
when the ^-component of the wave gets sufficiently large to 
reverse the flow near the surface. A lack of knowledge re- 
garding large-amplitude waves results in a lack of informa- 
tion concerning the point where transition will occur. 

It has not been proved that separation is a necessary 
condition for transition. Figure 11, where a slight reversal 
near the surface is indicated by a frequency doubling in the 
lower part of the cycle in the trace at 4 feet, gives evidence 
that reversal alone is not a sufficient condition. Perhaps 
the former picture, based on the analogy to the vortex sheet, 
is the more accurate one. 

Since the Taylor relation is not applicable when transition 
results from amplified waves, the average scale of the initial 
disturbance is important only insofar as it is related to the 
spectrum of the disturbance. The frequency distribution 
of the initial disturbance is of great importance, since 
amplification finally selects only a very narrow band of 
frequencies from all those initially present to bring about 
transition. The presence or absence of frequencies near this 
amplified band is obviously more important than the total 
energy of the disturbance. 

Only on the basis of certain prominent frequency bands 
in the stream turbulence can figure 7 be explained. The 
reasoning is as follows: If no disturbances were present, 
theoretically transition should not occur. The pairs of 
curves in figure 7 should then approach R x = 00 for zero 
turbulence. On the basis of amplification of initial dis- 
turbances in which all frequencies are present, equation (33) 
suggests that the curves defining the transition region should 
fall rapidly at first and then at a progressively diminishing 
rate as the turbulence is further increased. Actually, the 
opposite trend is shown in figure 7. Figure 8 shows that 
not all frequencies are represented alike and, furthermore, 
that noise is responsible for many of the peaks shown. As 
the present turbulence was decreased by the addition of 
damping screens, the ratio of noise to true turbulence was 
increased. While the total energy of the disturbance was 
decreased, the energy in certain frequency bands may actu- 
ally have been increased. A reduction of noise might have 
greatly altered figure 7. 


X. CONCLUSION 

Laminar-boundary-layer oscillations are the velocity 
fluctuations that result from a wave traveling downstream 
through the boundary layer. The characteristics of this 
wave have been measured and good agreement with theory 
has been obtained. The wave may be set up either by a 
vibrating object in the boundary layer, such as a vibrating 
ribbon, or by disturbances from the outside, such as stream 
turbulence and sound. When the disturbance is random or 
nearly random and contains many frequencies, amplification 
and damping isolate a wave containing a narrow band of 
frequencies in the neighborhood of the frequency most 
highly amplified. 

A wave in the boundary layer constitutes a disturbance 
that will cause transition on a flat plate with zero pressure 
gradient when the amplitude is sufficiently large. A suf- 
ficiently large random disturbance will also cause transition. 
When random initial disturbances are so small that transi- 
tion is delayed until values of R ^> 2000 have been reached, 
sufficient amplification can occur from the experimentally 
established lower limit 7?=450, to R 5 2000, to give rise 
to a recognizable oscillation. Larger initial disturbances 
obscure the oscillations by irregular fluctuations of com- 
parable magnitude. 

Natural disturbances may be sound waves as well as 
turbulence. T urbulence generally has a more nearly random 
distribution of energy with frequency than sound. In view 
of the importance of frequency, sound disturbances with a 
concentration of energy in frequency bands that are highly 
amplified by the boundary layer may be more conducive to 
early transition than turbulence. This is important in free 
flight where turbulence is probably negligible but where 
engine and propeller noise is present in large amounts. 

The effect of pressure gradient on the oscillations is a 
practicably important phase of the problem about which 
more should be known. Little more has been done here 
than to show the direction of the effect and to give some 
idea of its magnitude. The results are in accord with theory, 
which also is incomplete in this respect. 

Other aspects of the problem remain to be investigated. 
One of these having practical importance is the effect of 
curvature. The importance of boundary-layer oscillations 
on airfoils cannot be determined until effects of curvature as 
well as pressure gradient are known. Other quantities 
associated with the wave but not investigated are v and uv. 
A study of these quantities would be interesting as a further 
check on the theory, although hardly necessary to identify 
the oscillations and their characteristics in view of the 
information derived from studies of u. 

It is possible that boundary-layer oscillations may arise 
from internal disturbances as well as from external dis- 
turbances — that is, from surface irregularities and vibration 
of the surface. A randomly distributed small roughness 
may produce effects similar to small amounts of turbulence 
in the air stream. Small ridges or waves in the surface may 
start oscillations when the spacing is near some amplified 
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wave length. Vibration of the surface, like sound, may 
produce oscillations, especially when the frequency is near 
some amplified oscillation frequency. An investigation of 
these and other phases of the problem will throw additional 
light on the important problem of transition. 


National Bureau of Standards, 

Washington, D. C., April 1943. 
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Axis 

Force 
(parallel 
to axis) 
symbol 

Moment about axis 

Angle 

Velocities 

Designation 

Sym- 

bol 

Designation 

Sym- 

bol 

Positive 

direction 

Designa- 

tion 

Sym- 

bol 

Linear 
(compo- 
nent along 
axis) 

Angular 

Longitudinal— 

X 

X 

Rolling _ 

L 

Y >Z 

Roll 

0 

u 

V 

Lateral 

Y 

Y 

Pitching 

M 

z — >x 

Pitch. 

0 

v 

q 

Normal 

z 

z 

Yawing 

N 

X >Y 

Yaw 

xf/ 

w 

r 











Absolute coefficients of moment 

r _JL n 

Cl ~qbS Cm ~qcS 

(rolling) (pitching) 


Cn ~q,bS 

(yawing) 


Angle of set of control surface (relative to neutral 
position), 8. (Indicate surface by proper subscript.) 


4. PROPELLER SYMBOLS 


D 

P 

pID 

V' 

V. 

T 

Q 


Diameter 
Geotnetric pitch 
Pitch ratio 
Inflow velocity 
Slipstream velocity 

T 

Thrust, absolute coefficient C T = — tfu 

pn'Lr 

Torque, absolute coefficient <7 0 = — ^ 


C t 

v 

n 

$ 


Power, absolute coefficient 

V pV b 
Tn 2 

Efficiency 

Revolutions per second, rps 
Effective helix angle = tan - 


5. NUMERICAL RELATIONS 


1 hp =76.04 kg-m/s = 550 ft-lb/sec 
1 metric horsepower=0.9863 hp 
1 mph= 0.4470 mps 
1 mps=2.2369 mph 


1 lb=0.4536 kg 
1 kg=2.2046 lb 
1 mi= 1,609.35 m= 5,280 ft 
1 m=3.2808 ft 
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